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ABSTRACT 

We present optical spectra for a representative sample of 27 nearby (z < 0.2) 2MASS- 
selected AGN with red near-IR colours (J-K5 > 2.0). The spectra were taken with 
the ISIS spectrograph on the WHT with the aim of determining the nature of the red 
2MASS AGN, in particular whether they are young quasars obscured by their natal 
cocoon of gas and dust. We compare our findings with those obtained for comparison 
samples of PG quasars and unobscured type 1 AGN. 

The spectra show a remarkable variety, including moderately reddened type 1 
objects (45%), type 1 objects that appear similar to traditional UV/optical selected 
AGN (11%), narrow-hne Seyfert 1 AGN (15%), type 2 AGN (22%) and Hll/composite 
objects (7%). The high Balmer decrements that we measure in many of the type 1 
objects are consistent with their red J-K5 colours being due to moderate levels of dust 
extinction (0.2 < E(B-V) < 1.2). However, we measure only modest velocity shifts and 
widths for the broader [OIII]A5007 emission line components that are similar to those 
measured in the comparison samples. This suggests that the outflows in the red 2MASS 
objects are not unusual compared with those of optical/UV selected AGN of similar 
luminosity. In addition, the Eddington ratios for the 2MASS sample are relatively 
modest. 

Overall, based on their optical spectra, we find no clear evidence that the popu- 
lation of red, 2MASS selected AGN at low redshifts represents young quasars. Most 
plausibly, these objects are normal type 1 AGN that are moderately obscured by ma- 
terial in the outer layers of the circum- nuclear tori or in the disks of the host galaxies. 

Key words: galaxies: active — galaxies:Seyfert — quasars:emission lines — 
quasars:general 



1 INTRODUCTION 

Most of our understanding of Active Galactic Nuclei (AGN) 
is based on samples of AGN detected in UV/optical surveys. 
Indeed, key advances such as measuring high cosmological 



redshifts ( Hazard, Mackey & Shimmins 1963 1 and devel 



oping the orientation-based unified schemes (see Antonucci 
[T9931 ) are based on the study of AGN at UV/optical wave- 
lengths. However, there has been increasing evidence that a 
large fraction of AGN have not been detected in UV/optical 
surveys because they are obscured by large columns of dust 
( [Webster et al.|1995[[Cutri et al.|2002[ ). Surveys at infrared 
and radio wavelengths suggest that up to 80% of all AGN 
may have been missed in current UV/optical surveys (see 
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(MR); 



Low e t aI.'T988l [Webster et aL][T995l [Francis, Whiting fc| 



Webster„1999) . 



A key result of the Two Micron All Sky Survey 
(2MASS) was the discovery of a population of AGN that ap- 
pear redder than their traditional optical/UV selected coun- 
terparts at infrared wavelengths. The red 2MASS AGN were 
selected to have J-Kg > 2.0, distinguishing them from most 
optical/UV selected AGN that have bluer near-lR colours 
(J-Ks < 2.0; [Cutri et aL][2002[), and K-band magnitudes 



in the range 11.0 < Ks < 14.9 mag. The latter criterion 
was chosen so that bright (Kg < 11.0) Galactic AGB and 
carbon stars were excluded, along with 'normal' galaxies in 
the redshift range 0.4 < z < 0.5 (Kg > 14.9; jCutri et al.[ 
2002). Spectroscopic follow up of the objects which fulfill 
these criteria revealed that the majority of the candidates 
have spectra resembling AGN ( [Cutri et al.[[2002| ). More re- 
cent surveys have determined that these 'red quasars' do 
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indeed represent a large fraction of the overall population 
of quasars in the local Universe (~15-60%; Glikman et al 
2007| [Glikman at al.||2012[ ) 



list of Hutchings et al. (20031, which is itself representative 



of the population of red, 2MASS-selected quasars. We also 



Because the red quasars potentially make up a large 
fraction of the AGN population, it is important that we un- 
derstand their relationship to the unobscured population: 
are they young, dust obscured versions of 'typical quasars' 
([Hu tchings ct al. 2003; Gcorgak akis et al.f 2009': 'Glikman eti 
|al.|2007,,Glikman at al.,2012, )? Or can their properties be ex- 
plained by a specific viewing angle in the orientation-based 



unified schemes (Wilkes et al. 20021? Or do they perhaps 



represent a new population altogether? Are these objects all 
luminous quasars, or are some less luminous AGN? 

Past studies of red 2MASS quasars have found redden- 
ings in the range 0.1 < E(B-V) < 3.2 magnitudes based on 
Balmer decrements and fits to the continuum shape (e.g. 
Glikman ct al. 2 007{ [Georgakakis et al..2009: Kuraszkiewicz| 



et al.||2009a,b] [Glikman at al.||2012[ |Canalizo et al.||2012[ )" 



although it is important to recognise that some have negh- 
gible reddening ( Glikman et al.|2007 |. In addition, imaging 
studies of the host galaxies find evidence for disturbed mor- 
phologies in up to 70% of cases (Hutchings et al.|2003 Mar 



ble et al. 2003 Urrutia, Lacy & Becker 2008), which may 



support the idea that host galaxies have undergone recent in- 
teractions with other galaxies. Such interactions could drive 
material into the central region of the galaxies, enshroud- 
ing the SMBH in a cocoons of dust (Hopkins et al. 2006) , 
making them difficult to detect in traditional UV/optical 
surveys ( [Cutri et al"[[2002[ ). 

Many previous studies of 'red quasars' have focussed 
on samples covering a wide range of redshift (0 < z < 3: 
Glikman et~aL[[2007[ [Georgakakis et al.|[2009i [Glikman at"[ 



al.|2012 1. The problem with such samples is it is difficult to 



investigate the intrinsic diversity of the red quasar popula- 
tion in the face of possible evolutionary and/or luminosity- 
dependent effects. Therefore it is important to study a sam- 
ple of 2MASS-selected AGN that covers a more limited range 
in redshift. 

This paper reports observations of an RA-limited sam- 
ple of 27 2MASS selected quasars (J-Kg > 2.0) with a nar- 
row range of redshifts: 0.09 < z < 0.20. We aim to deter- 
mine how the red 2MASS AGN fit into our picture of active 
galaxies in general, and how important they are to our un- 
derstanding of AGN. In this, the first of two papers (Paper 
I), we report a spectroscopic study of this sample, focussing 
on the AGN properties, emission line outflows and dust ex- 
tinction. The spectral properties are compared to control 
samples of optical/UV selected quasars and AGN. Paper II 
in the series will focus on the near- and mid-IR properties 
of the sample, and will include a comprehensive comparison 
with several complete samples of active galaxies. 

The cosmological parameters used throughout this pa- 
per are adopted from WMAP: Hp = 71 km s"\ Am = 0.27 
and = 0.73 ( [Spergel et al][2053 l. 



2 OBSERVATIONS AND DATA REDUCTION 

The full sample for this survey comprises a complete RA- 
limited (2 < RA < 15 hr) sub-sample of 27 objects with 
red near-IR colours (J-K5 > 2), K-band magnitudes 11.0 
< Ks < 14.9 mag, and redshifts z < 0.2 selected from the 



observed a further 2 objects from Hutchings et al. ( 2003 1 
that fall outside the original redshift and /or RA range: 
J1637+25 (z=0.277), and J2124-17 (z^O.lll). Therefore, 
the total sample consists of 29 objects, however two of these 
objects have redshifts z>0.2. Although we present the data 
for the full 29 objects, we only consider the results for the 27 
objects with redshifts z<0.2 in the analysis. Table [l] shows 
the basic properties of the sample. 

Low resolution optical spectroscopic observations of the 
sample were taken with the ISIS dual-arm spectrograph on 
the 4.2-m Wilfiam Herschel Telescope (WHT) on La Palma 
in three runs: the 28th of July 2006, 8th and 9th of Febru- 
ary 2007, and the 26th and 27th of September 2011. The 
data were taken with the ISIS Dual arm spectrograph on 
the WHT. On the red arm, the R158R grating was used 
with the REDPLUS CCD, and on the blue arm, the R300B 
grating was used with the EEV12 CCD. A dichroic cutting 
at 53OOA was employed to obtain spectra simultaneously 
in the blue (~3250-5250A), and in the red (~5200-9500A). 
To reduce the effects of differential refraction, all exposures 
were taken when the objects were at low airmass (sec z < 
1.1) and/or with the slit aligned close to the parallactic an- 
gle. The seeing for the nights of the observations varied over 
the range 0.8 < FWHM < 1.3 arcseconds. 

Sets of three exposures were taken on both arms si- 
multaneously for each object using a 1.5 arcsecond slit. An 
additional wider slit (5 arcseconds) exposure was taken on 
both arms for 23 of the object^ This was to assess the effect 
of possible slit losses on the emission line flux. To eliminate 
contamination from second order emission, a GG495 block- 
ing filter was introduced into the ISIS red arm. The exposure 
times for each object are given in Table [T] 

The data were reduced in the standard way (bias sub- 
traction, fiat fielding, cosmic ray removal, wavelength cali- 
bration, flux calibration) using packages in iRAlj^ The two- 
dimensional spectra were also corrected for spatial distor- 
tions. To reduce wavelength calibration errors due to flex- 
ure of the telescope and instrument, arc spectra were taken 
at the position of the object on the sky. Based on mea- 
surements of the night sky lines, the wavelength calibration 
uncertainties are 0.1 A on both the blue and red arms, and 
the narrow slit spectral resolution is estimated to be ~5.6A 
and ~10.5A on the blue and red arms respectively for 23/29 
of the objects, and ~6.6A and ~6.lA for the blue and red 
arms for the rest of the sample (see Table [l|. Atmospheric 
absorption features were removed by dividing the red spec- 
trum of each object by that of a telluric standard, taken 
close in time and airmass to the observations of 27/29 of 
the objects. No telluric standards were taken for J1637-I-25 
and ,12124-17, and therefore the atmospheric features were 



^ Wide slit exposures were not taken for 6 of the objects in 
the sample: J0221+13, J0248-I-14, J0306-05, J0312-I-07, J1637+25 
and J2124-17. These objects were observed on separate runs to 
the rest of the sample and no wide slit exposures were taken for 
them. 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
the Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation ( |http:/ /iraf.noao.edu/| . 
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Figure 1. The combined (blue and red arm) spectra of tlie 2MASS sample discussed in this paper. Note the remarkable variety in the 
spectra: the spectral types range from Seyfert 2 galaxies that lack the broad emission lines (e.g. J0504-19 & J1057-13), through highly 
reddened broad line objects (e.g. J1040+59 & J1127+24), to objects that show optical spectra similar to 'normal' blue quasars (e.g. 
J1338-0438 & J2124-17), peculiar objects with strong forbidden high ionization lines (FHILs) in their spectra (e.g. J1131+16; |Rose et al.| 
|2011[ | and even LINER/HII region-like objects (e.g. J0411-01 & JlOOl-l-41). J1637+25 and J2124-17 have not had the telluric absorption 
features removed, because no telluric standard stars were observed during the particular run in which they were measured. Continued 
ON NEXT PAGE. 



not removed for those objects. The spatial pixel scales of the 
2D spectra are 0.4 arcseconds in the blue and 0.44 arcsec- 
onds in the red, and the relative flux calibration uncertainty 
- based on 11 observations of 8 flux standard stars taken 
throughout the runs - is estimated to be ±5%. 

The spectra were extracted and analysed using the 
STARLINK packages FIGARO and dipso. All extracted spec- 
tra were corrected for Galactic extinction using the Galac- 



tic extinction maps of Schlegel, Finkbeiner & Davis ( 1998 1 



and the extinction laws of of Seaton ( 1979 1 for the UV, and 
|Howarth| ( |1983[ ) for the optical/IR, prior to the analysis. In 
order to investigate the nuclear spectra, we extracted aper- 



tures centred on the nuclei of the targets in the 2D spectra, 
with aperture sizes chosen such that the full flux from the 
nucleus of each object was extracted. The aperture sizes are 
given in Table [T] and the extracted spectra are shown in 
Figure [l] 
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3 SAMPLES AND EMISSION LINE FITTING 



3.1 2MASS sample 

The basic properties of the 2MASS sample are given in Ta- 
ble [l] To give an indication of the intrinsic powers of the 
AGN, the [OIII] emission line luminosities ijjyoiii]) were 
calculated using fluxes measured from the wide slit data for 
most objects; the narrow slit data were used where wide 
slit data was not available. For the purpose of calculating 
L[o///l the [OIII] fluxes were corrected for Galactic redden- 
ing but not intrinsic reddening caused by dust in the host 
galaxies, hyoiii] was chosen because the [OIIIJA5007 emis- 
sion line is regarded as a good indicator of AGN bolometric 



powe r ( |Reyes et al.|2008|[D"icken et al.|2009[|LaMassa et al 

20101, and can be used to distinguish between quasars and 
Seyfert galaxies (IZakamska et al.|2003[ ). The redshifts of the 
objects were determined by fitting single Gaussian profiles 
to the [OIII]A5007 emission lines. 

The wide-to-narrow slit [OIII] flux ratios are shown in 
Figure [2] The majority of objects have wide and narrow slit 
fluxes that agree within ±40%, suggesting that the narrower 
slits captured the majority of the flux from the emission re- 
gions: the median wide to narrow flux ratio is 1.14, showing 
that the wide slit captured slightly more flux than the nar- 
row slit, as expected because the wide slit data probes a 
larger physical scale (4.9 to 22.7 kpc for the redshift range 
of this data) when compared to the the narrow slit data (1.5 
to 6.8 kpc). There are two outliers: J0435-06 and JlOOl-l-41, 
which have ratios above 1.5. However, the wide to narrow 
flux ratios for these objects agree within ~2cr of the median 
ratio. In these cases the [OIII]A5007 emission lines have a 
low equivalent width. Therefore the large wide to narrow 
flux ratios in these objects may reflect the large uncertain- 
ties in [OIII] fluxes, as measured from the noisy wide slit 
data, rather than uncertainties in the narrow slit fluxes. The 
similarity between the flux captured in the wide and narrow 
slit observations gives confidence that the spectra for the 
5 objects without wide slit data (runs 'b' and 'c' in Table 
pi give an accurate indication (within ~40%) of L[oi//] for 
these objects. 



3.2 Comparison Samples 

In attempting to understand the nature of the red 2MASS 
AGN, some of the key questions we are trying to address 
are: if one takes a sample of AGN selected on the basis of 
their red colours at near-IR wavelengths (J-K>2.0), how do 
the optical spectroscopic properties of such a sample dif- 
fer from those of samples of typical type 1 AGN selected 
using their optical/UV colours, and can any differences be- 
tween the samples be explained in terms of reddening ef- 
fects? Therefore, we have compared our results with two 
comparison samples of 'typical' quasars. For any compar- 
isons we make, the samples we use are limited to redshifts 
z<0.2, so they are consistent with our 2MASS AGN sample. 



3.2.1 PG Quasars 

One of the comparison samples used throughout this paper 
is the PG quasar sample of |Boroson fc Green] ( [19921 ) . The 
PG quasars comprise a complete sample of 87 UV-selected 
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^1.5- -r 



0.5 




Figure 2. Wide to narrow [OIII] flux ratio plotted against red- 
shift for the 2MASS sample. 



objects from the Palomar bright quasar survey (Schmidt 
& Green 19831, selected to have redshifts z < 0.5, and 



UV/optical colours U-B < -0.44 to ensure an excess of fiux 
in the UV. This sample is also selected to have quasar-like 
luminosities, with a total brightness of Ms < -23 for the 
AGN and host galaxy together, and broad permitted emis- 



sion lines clearly present in the optical spectra (Boroson 
fc Green| [T992). The PG quasars provide a useful compari- 



son with the 2MASS sample because they represent 'typical' 
blue selected quasars at their redshifts (z < 0.5). By limiting 
the redshift range of the sample, there are 61 PG quasars 
suitable for comparison with the red 2MASS sample. Spectra 



of the Boroson & Green ( 1992 1 sample were kindly supplied 



by T. Boroson, and remeasured by MR. The data cover a 
spectral range of AA4300-5700 and have spectral resolutions 
in the range 6.5-7.0 A. Detailed properties of the sample 
and the Fell subtraction procedure are outlined in |Boroson| 
|fc Green] ( p92l ). 

The spectra of the PG quasars were fitted in a similar 
manner to the 2MASS objects (see ^ 3.3). However, in con- 
trast to the 2MASS sample, only H/3 and [OIII]AA5007,4959 
were fitted for the PG quasars, due to the limited spectral 
range of their spectra ( |Boroson Sz Green|1992 |. 

A sub-sample of these quasars (17) have high quality 
Sloan Digital Sky Survey (SDSS; ]York et al.]|2000l ) spectra 
available, which cover a wider wavelength range than the 
spectra from the ]Boroson fc Green ( 1992 1 study. We will 
refer to this as the PG/SDSS sample. Where possible, these 
spectra have been used as part of the comparison with the 
2MASS sample. 

As shown in Figure ]3j which plots the J magnitude 
against the J-Kg colours, the PG quasars are compara- 
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The spectral fitting teclinique outlined in |Jin et aL] 



2.5 




1.5 - 



Figure 3. J-Kg versus J for the 2MASS and comparison samples. 
The 2MASS population measurements are indicated by the filled 
circles, the PG quasar population l |Boroson &c Green||1992] | are 
indicated by the unfilled squares and the unobscured type 1 AGN 
( |Jin et al.|2012a,b,c[ l are indicated by asterisks. Higher values of 
J-Ks indicate redder colours. 



tively 'bluer' than the 2MASS sample at near-IR wave- 
lengths (PG median: J-Ks = 1.79±0.03; 2MASS median: J- 
Ks=2.15±0.08), with minimal overlap in J-Kg colours be- 
tween the samples: only 20% of the PG sample have J-Ks 
> 2.0. In addition, the 2MASS objects are fainter in the J 
band. Given that the redshifts of the 2MASS objects in our 
sample are lower on average than those of the PG quasars, 
this suggests that the 2MASS objects have lower luminosi- 
ties and therefore we may not be able to consider them as 
true quasars (see f 4.1.1 for a full discussion). 



3.2.2 Unobscured Type 1 AGN 

The second comparison sample comprises a sample of 51 
relatively nearby (z < 0.38), unobscured type 1 AGN se- 
lected by |Jin et aT] ( |2012a,b^ on the basis of low optical 
reddening, and low gas columns as indicated by their X-ray 
spectra. The optical properties of this sample were obtained 
from high quality SDSS DR7 spectra ( |Jin et al.|2012a,b,c[ ). 
The SDSS DR7 spectra cover a wavelength range of 3800- 
9200A and have spectral resolution of ~3A. The analysis 
techniques are outlined in the work of |Jin et al.| ( |2012a,b,c[ ). 
There are 14 PG quasars in this sample, therefore we have 
removed these from any comparison which includes the PG 
quasar sample to avoid overlap. In addition, when consider- 
ing the remaining objects with redshifts z<0.2, there are 21 
unobscured type 1 AGN suitable for comparison. As these 
objects comprise unobscured AGN, we consider this a suit- 
able sample to compare to the red 2MASS quasars. 



(2012a I is similar to the methods we have used to fit the 
emission lines of our 2MASS sample (see ^ 3.3). Therefore 
we are confident that we can use the data presented in | Jin| 



et al. (2012a I for a direct comparison with our results. 

As shown in Figure |3j like the PG quasars, the un- 
obscured type 1 AGN are also relatively 'blue' compared 
with the 2MASS sample at near-IR wavelengths (median J- 
Ks=1.60±0.04), with no overlap in J-Ks colours between 
the samples. Unlike the PG quasars, this sample has com- 
parable J-band magnitudes to the 2MASS sample. This sug- 
gests that, given the redshifts of the objects in the sample 
(median z~0.16; for the entire sample), this sample has sim- 
ilar luminosities to the 2MASS objects. 



3.3 Continuum Subtraction and Emission Line 
Modelling 

Fitting emission features in AGN spectra is often compli- 
cated by the presence of the host galaxy continuum in the 
spectrum. Subtle features such as weak broad components 
can often be overlooked, and the presence of underlying 
Balmer absorption from the stellar continuum can lead to 
underestimation of the intensity of the Balmer and other 
emission line strengths (e.g. Sarzi et al.|2006| ). 

For the purposes of studying the emission regions, the 
underlying continuum was modelled and subtracted in all 
objects in which the stellar continuum is significant — 
mainly LINER/Seyfert 2 type objects (e.g. J0221-I-13) and 
objects with weak broad lines (e.g. J0422-18). These steps 



required the use of the starlink package dipso ( Howarth et 



called CONFIX ( 


Robinson et al.|2000 


Tadhunter et al., 20051 


Holt et al.||2007| |Rodn'guez Zaurfn, Tadhunter & Gonzalez] 


Delgado|2009 


I. 



We started by experimenting with subtracting a neb- 
ular continuum component, which might be significant at 



near-UV wavelengths (Dickson et al. 19951. The contribu- 



tion of the nebular continuum was estimated using the nar- 
row H/3 emission line flux. However, its strength is sensi- 
tive to the reddening of the NLR, which is not accurately 
known in most cases. Indeed, we found that in all cases we 
obtained better fits if we did not subtract a nebular contin- 
uum. Therefore the results below refer to models without 
nebular continuum subtraction. 

Once the spectra had been prepared, the continua 
were modelled using the customised idl procedure confix 
( Robinson et al.|2000||Holt et al.|2007[ ). To minimise degen- 
eracy when fitting our spectra, we fit the smallest number 
of stellar and/or power-law components to the continuum as 
possible. For each object we ran 3 sets of models. In the first 
run we fit an unreddened old stellar population (OSP, 12.5 
Gyr) and a young stellar population (YSP), with ages in 
the range 0.001-5 Gyr, and reddening: < E(B-V) < 2.(0 
The inclusion of a YSP is justified on the basis that many 
of the objects with significant stellar continua show higher 



3 We used the |Calzetti et al.| | [2000l l reddening law to redden the 
YSP spectra. 
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Figure 4. A best-fit CONFIT model for J10404364-I-5934092 (x^=0.54), plots (a)-(d) show how the combined stellar populations fit the 
detailed features. The object requires an OSP, a YSP, and a power-law component for an adequate fit. The YSP has an age of 0.005 
Gyr, an extinction of E(B-V)=1.0 and represents 58.7% of the flux in the normalising bin. The power-law component has a=5.11 and 
represents 16.4% of the flux in the normalising bin. Finally, the OSP represents 24.8% of the flux in the normalising bin. The flux is in 
units of erg s~^ cm~^ A"-*^, and wavelength is in units of A. 
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Table 2. The [OIII] model parameters used when fitting the emis- 
sion lines of the spectra. The second and third columns present 
the rest-frame widths of the narrow and broad kinematic compo- 
nents respectively. The final column gives the rest-frame velocity 
separation of the broad relative to the narrow kinematic compo- 
nents. 



Object 


FWHMjv 


FWHMs 


Separation 




(km s^-"-) 


(km s^-"-) 


(km s^^) 


J0221-I-13 


Inst. 


799 it 36 


-217 it 27 


J0248-I-14 


448 it 39 


945 it 61 


-168 it 44 


J0306-05 


Inst. 


826 it 148 


-426 it 152 


J0312-f07 


Inst. 


718 it 95 


-227 it 72 


J0400-f05 


326 it 13 


1197 it 51 


-36 it 11 


J0409-I-07 


620 it 22 


1842 it 340 


-44 it 68 


J0411-011° 


544 ± 150 






J0422-18'' 


Inst. 


1700 it 56 


-298 it 31 


J0435-06 


256 it 19 


376 it 84 


-655 it 93 


J0447-16 


194 it 7 


1072 it 165 


-569 it 187 


J0504-19'' 


Inst. 


1778 it 40 


191 it 32 


J0910-I-33 


715 it 21 


1652 it 48 


-276 it 33 


JlOOH-41'' 


418 it 76 






J1006-f4l'' 


Inst. 


554 it 159 


-353 it 152 


J1014-I-19 


428 it 13 


1637 it 63 


127 it 23 


J1040-I-59 


181 it 22 


993 it 84 


214 it 62 


J1057-13 


707 it 14 


1477 it 139 


-855 it 116 


J1127-I-24 


267 it 7 






JllSl-l-ie^ 


Inst. 


720 it 31 


-92 it 5 


J1158-30 


693 it 27 


1177 it 139 


-1195 it 66 


J1212-14'' 


455 it 25 






J1321-I-13" 


390 it 14 






J1323-02 


366 it 16 


1512 it 108 


-488 it 94 


J1338-04'' 


Inst. 


780 it 24 


-224 it 34 


J1407-f42 


584 it 8 


1316 it 48 


190 it 28 


J1448-I-44" 


441 it 15 






J2124-17" 


753 it 22 







" The model for this object requires just one component. 
^ The narrow component of the model can be modelled with 
the instrumental width. 
This object was studied in detail in |Rose et al.|p011| . 

order Balmer lines in absorption (e.g. JlOOH-41). The stel- 
lar population templates were taken from the instantaneous 



Chariot (20031 



burst spectral synthesis results of Bruzal 
for a Salpeter IMF and solar metallicities. For the second 
run a power-law component, which represents either the di- 
rect or scattered AGN component, was also included in the 
model fitting. Finally, an OSP with and without a power-law 
was fitted to each spectrum, independent of a YSP. 



The spectra were fitted using fluxes measured in 30A 
bins, chosen to have as even a distribution in wavelength 
as possible, and avoiding the brighter atmospheric emission 
and stronger atmospheric absorption, and AGN emission, 
features in the spectra. The entire useful wavelength ranges 
of the spectra were fitted: observed wavelengths 3200- 
9000A. In addition, a normalising bin was chosen which is 
also free of any emission/absorption features (~5050-5200 
A for all objects in the rest frame). 

Fits to the overall continuum s hape were chosen with 
low reduced <1 (see discussion in 



Tadhunter et al. 



20051 



Then a visual inspection of fits to individual absorption 



lines, including Call K (3934A), Mgl (5175 A), Balmer and 
Nal D (5890-5896 A), was used to select the best overall fits. 
Note that not all objects were successfully fitted. For exam- 
ple, we could not produce an adequate fit to the type 2 object 
J113H-16 because there are relatively few genuine contin- 
uum bins suitable to model the spectrum of this object. It 
should be noted that this object has an unusual spectrum, 
which contains a large number of forbidden high ionisation 
lines (FHILs; see Rose et al. 2011[ |. Figure [4] shows an ex- 
ample of a successful confit fit to a 2MASS object. Details 
of the models chosen for the continuum subtraction in indi- 
vidual objects are given in Table ID 

Once the continuum had been subtracted from the spec- 
tra (where applicable), the emission lines were fitted using 
the STARLINK DIPSO package (Howarth ct al. 2004). 

The emission lines in the spectra of the 2MASS and 
PC comparison sample AGN were initially fitted with single 
Gaussian profiles, starting with [OIII]AA5007,4959 doublet. 
We used the [OIII] AA5007, 4959 emission lines because they 
are often the strongest unblended lines emitted by the NLR. 
To produce the models, Gaussian profiles were fitted to the 
[OIII] doublet, with the separation and intensity ratio set by 
atomic physics for each kinematic component (as defined by 
the velocity shifts and line widths). While single Gaussian 
fits were sufficient in some objects, such fits did not ade- 
quately model the wings of the emission lines of most of the 
objects; in such objects double Gaussian model were fitted to 
the spectral lines. When a second Gaussian was added, the 
line ratios, wavelength separations and FWHM were fixed 
for the two doublet components as before, however the shift 
and FWHM were allowed to vary relative to the first compo- 
nent, and the widths were allowed to vary. In all cases, one 
or two Gaussian components proved sufficient, and in all 
cases where two Gaussian profiles were required, the width 
of one Gaussian was significantly broader than the other. 
In some cases, the widths of the narrow components for the 
[OIII] emission lines were found to be consistent with the 
instrumental width of their spectra. In such cases, the nar- 
row components of the models that were fitted to the other 
lines were fixed to the instrumental width. However, in the 
other (resolved) cases, the measured FWHM was used to 
obtain an intrinsic velocity width for the narrow component 
by correcting its FWHM in quadrature using the instrumen- 
tal width. The same technique was applied to the broader 
components of the model. Also, the broad component is of- 
ten shifted with respect to the narrow component; this is 
accounted for in the models. 

Once the model parameters had been determined from 
the fits to [OIII], they were used to fit the other lines 
in the spectra. The [OIII] model parameters for each ob- 
ject are presented in Table [5] Known blends (e.g. Ha & 
[NII]AA6548,6583) and doublets (e.g. [SII] AA6716,6731 and 
[OII]AA3726,3729) were modelled with the [OIII] model 
using constraints provided by atomic physics (i.e. fixed 
FWHM, line separation and, where appropriate, the inten- 
sity ratios set by atomic physics). In most cases, the [OIII] 
model successfully fit all the emission lines/blends in the 
spectra. However, in a few cases where emission lines have 
a low equivalent width (e.g. weaker lines such as [OI] or 
[FeX]), one of the components in the 2 component [OIII] 
model provides a better fit to the emission line. This hap- 
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Table 3. Best fit CONFIT model results for the 2MASS sample. The simple stellar population models were taken from |Bru2al fc Charlot| 
l |2003[ |. The subscripts 'U' and 'L' refer to the upper and lower cases of the model parameters that provide a good fit to the continuum. 
For the cases where the best fits include power laws we use the form oc v'^°' , where a represents the power law index. 



Object YSPlt (Gyr) YSPl (Gyr) E{B-V)y E{B-V)£ au Ql 



J0221+13 


1.0 


0.7 


0.5 


0.3 


2.18 


0.457 


J0306-05 


2.0 


0.8 


0.5 


0.3 


3.95 


1.40 


J0411-01 


0.006 


0.005 


1.1 


0.6 






J0422-18 


0.005 


0.004 


1.1 


0.9 


2.35 


0.32 


J0447-16 


1.2 


0.8 


0.3 


0.2 


1.31 


1.28 


J0504-19 


0.01 


0.006 


0.9 


0.7 






J0910+33 


0.9 


0.3 


0.3 


0.1 


-0.347 


-0.818 


JlOOl+41 


0.2 


0.1 


0.7 


0.6 






J1040+59 


0.005 


0.004 


1.1 


0.8 


5.97 


5.11 


J1057-13 


2.0 


0.2 


0.6 


0.0 


-0.30 


-1.33 


J1158-30 


0.4 


0.2 


0.5 


0.2 


2.86 


1.74 


J1212-14 


1.2 


0.8 


0.4 


0.2 


1.07 


0.627 


J1321+13 


0.05 


0.005 


0.7 


0.4 






J1407+42 


0.06 


0.04 


0.6 


0.5 







pened in some objects with strong broad Balmer emission 
lines. 

The broad permitted lines emitted by many of the ob- 
jects (mainly Ha and H/3) were initially fitted with a sin- 
gle Gaussian component, where the intrinsic velocity widths 
(FWHM) were kept the same for the all BLR emission lines 
in each object. The NLR emission components of the Balmer 
emission were still fitted with the [OIII] model. However, for 
11 objects in the 2MASS sample, the broad and narrow line 
components of the Balmer lines could not be confidently sep- 
arated because, at the resolution of our spectra, they merge 
seamlessly together, leading to degeneracies in fits that at- 
tempt to model the lines as single narrow and broad Gaus- 
sians (see Figure[5|. In addition, for the PG/SDSS sample, 
we could only confidently separate the broad and narrow 
emission in 6 objects. Figure [5] highlights the problem of de- 
generacy in the blending of the NLR with BLR in the Balmer 
lines, using the example of J0422-18. Although the overall 
fit to the E[/3 feature is adequate, the narrow line Gaussian 
fitted to the H/? emission line may be over- or underestimat- 
ing the flux of the NLR If/3 emission. In such cases, we only 
consider the total fluxes of the combined broad and narrow 
emission components of the Balmer lines. 

The H/3 and [OIII] AA5007,4959 may also be strongly 
blended with broad permitted Fell emission. In cases where 
broad permitted Fell emission is prominent in the spec- 
tra (e.g. J0248-I-14), the Fell emission was fitted separately 
from the affected emission lines as follows: the Fell emission 
lines in the same multiplets (F, S or G) were fitted with 
the same intrinsic velocity widths as the broad H/3 emission 
line, and their expected intensity ratios were used to con- 
strain the fits, as outlined in'Kovacevic, P opovic fc Dimitri-| 
|jevic| ( |2010[ ). This approach provided a good fit to the Fell 
emission blends in all objects where they are significant. 

Examples of model fits are given in Figures |6] and [T] Fig- 
ure [6] shows an example of a multi-component Gaussian fit 
to the R/3 and [OIII]AA5007,4959, where the broad and nar- 
row components of H/3 were easy to separate. In contrast. 
Figure [7] shows a case where the broad and narrow emis- 
sion could not be easily separated, and where Fell emission 



Q0422-1854 - Aperture 1 




Wavelength 



Figure 5. Gaussian fits to the RI3 and [OIII] AA5007,4959 emis- 
sion lines of J0422-18. The overall fit is represented by a red line. 
Both the NLR [OIII] model (orange lines; the narrow component 
is represented by a dashed line and the broad component is pre- 
sented by a dotted line) and BLR components (represented by a 
blue dashed and dotted line) are displayed. The NLR model fit 
to the H/3 blend highlights the degeneracy issues with blended 
emission from the BLR and the NLR (this is described in more 
detail in ^ 3.3). 

is prominent. Interestingly, the BLR of this object is sub- 
stantially blueshifted relative to the NLR in this object (by 
-2870±160 km s~^). The properties of the individual objects 
in the 2MASS sample with particularly interesting spectra 
are discussed in the Appendix. 

The results for the NLR [OIII] model fits to the other 
lines are presented in Table [4] all the emission line ratios are 
presented relative to [OIII]A5007, and the narrow and broad 
line region Balmer components are presented separately. 
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Figure 6. Gaussian fits to ttie and [OIII] AA5007,4959 emis- 
sion lines of J0400+05. Both tlie NLR [OIII] model and BLR 
components are displayed. In this case the BLR and NLR com- 
ponents of H/3 could be clearly separated. Interestingly, the BLR 
is substantially blueshifted relative to the NLR in this object (by 
-2870itl60 km s~^). The properties of the individual objects in 
the 2MASS sample with particularly interesting spectra, are dis- 
cussed in the Appendix. 



QS02124 - Aperture 1 




Wavelength 



Figure 7. Gaussian fits to the and [OIII] AA5007,4959 emis- 
sion lines of J2121-17. The NLR [OIII] model, BLR component 
and Fell emission (represented by a violet dashed and dotted line) 
are displayed. The overall fit is far better than when the emission 
lines are modelled without the Fell emission. 
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tinuum and host galaxy contamination. Following |Zakamska| 
et al. (20031 we assume that an object is a quasar if it has 




Figure 8. [OIII] emission line luminosity {^\oill]) plotted 
against redshift for the 2MASS objects and the comparison sam- 
ples. The 2MASS objects are represented by filled circles, the PG 
quasars are represented by unfilled squares, and the unobscured 
type 1 AGN are represented by asterisks. Objects with l^yoiii] 
greater than 10^^ W are likely to have quasar-like bolometric lu- 
minosities. 



4 RESULTS 

The optical spectra of the 2MASS sample, shown in Figure 
[l] reveal a remarkable variety, ranging from Seyfert 2 galax- 
ies that lack broad emission lines (e.g. J0504-19 & J1057-13), 
through highly reddened broad line objects (e.g. J1040-(-59 
& J1127-I-24), to objects that show optical spectra similar to 
'normal' blue quasars (e.g. J1338-04 & J2124-17) and even 
LINER/HII region-like objects (e.g. J0411-01 & J1001-H41). 
We now describe the results obtained from the detailed anal- 
ysis of the spectra. 



an [OIII]A5007 luminosity Ijyoiii] > 10 W, corresponding 
to a nuclear continuum luminosity of Mb < -23 ( [Schmidt fc| 
Green||1983[ i Zakamska et al.|2003[ ). 

Figure 8| shows hyoiii] plotted against redshift for the 
2MASS objects and the comparison samples. Note that the 
luminosities have been determined using the observed frame 
[OIII]A5007 fluxes and the luminosity distances appropriate 
for the assumed cosmology, but have not been corrected for 
intrinsic reddening. Where possible we use the wide slit (5") 
data to ensure that the AGN [OUT] luminosity measure- 
ments are not affected by slit losses; for the 5 objects for 
which there are no wide slit data, we use the available 1.5 
arcsecond data. We find that 5 objects in our sample (19%) 
qualify as quasars according to the [OIII] emission line lumi- 
nosity criterion: J0400-I-05, J1131-I-16, J1338-04, J1407-f42 
and J2124-17. Note that the latter object was previously 



identified as a ULIRG by Kim fc Sanders ( 1998 1 



The rate of quasars for the 2MASS sample (filled circles) 
is comparable to that of the rate of quasars in the compari- 
son samples: 28% of the PG quasars (unfilled squares) and 
10% of the unobscured type 1 AGN (asterisks) have quasar- 
like luminosities based on the Ljo///] luminosity. Therefore 
the rate of occurrence of quasars in the PG and unobscured 
type I samples is not significantly different from that of 
the 2MASS sample. However, we note that 4 objects in our 
2MASS sample have lower [GUI] emission line luminosities 
than any of the objects in the PG quasar and unobscured 
type 1 AGN samples. 

We can only be confident that 5 of the objects in 
the sample are genuine quasars. The remaining sample has 
[OIII] luminosities in the range 32.5 < \^[oiii] < 35 W — 
more typical of Seyfert galaxies. If the 4 objects with low 
[OIII] luminosities are ignored, the 2MASS AGN have a 
similar range of luminosities as the PG quasars and the un- 
obscured type 1 AGN. Note that we have not corrected the 
[OIII] emission line luminosities for intrinsic dust extinc- 
tion. If we correct the emission line luminosities using the 
NLR reddening estimates available for 16 objects (^ 4.2.1), 
we find that the rate of occurrence of quasars in the 2MASS 
sample increases to 27% (8 objects). However, we emphasise 
that we can only accurately correct the narrow line fluxes 
for reddening in 55% of the objects. 



4.1 Quasars, AGN or starbursts? 

4.1.1 Quasars? 

Previous studies of red 2MASS AGN have shown that 
most red 2MASS AGN in the local Universe are not 
quasars, but have luminosities more consistent with those 
of Seyfert galaxies dSmith et al.||2002| |Marble et al.||2003 



Kuraszkiewicz et al.l2009a,b I. 



To test the incidence of quasars in our low-z 2MASS 
sample, we used the [OIII]A5007 emission line luminosity, 
which has been shown to be a good indicator of AGN bolo- 
metric power (see [Heckman et al. [2004 , Bian et al.|[2006[ 



Dicken et al.|[2009 LaMassa et al 



2010 1. Other luminosity 



indicators, such as the L5100 optical continuum luminosity, 
cannot be accurately calculated in most of the 2MASS ob- 
jects because of a combination of reddening of the AGN con- 



4.1.2 AGN or Star-bursts? 

Both starbursts and AGN activity can potentially ionise the 
narrow emission regions in the nuclei of galaxies ( [Baldwin^ 
Phillips & Terlevich 19811. It is therefore important to es- 



tablish whether the NLR in the 2MASS objects are indeed 
photoionized by AGN, or whether recent starburst activity 
provides the photons needed to ionize the emission line re- 
gions. Therefore we have categorised the spectra based on 
the emission line ratios using Baldwin, Phillips & Terlevich 



(BPT) diagnostic diagrams (e.g. see|Baldwin, Phillips & Ter- 
levich|198l||Kewley et al.||2006[) 



Figures |9|- |ll|show the results for the 2MASS sam- 
ple. The divisions between the different classes of objects in 
these figures are taken from Kewley et al. (20061. Note that 
only ~50% of the 2MASS sample could be plotted on the 
BPT diagrams because of the difficulty of measuring accu- 



16 M. Rose, C. N. Tadhunter, J. Holt, J. Rodriguez Zaurin 




Figure 9. Diagnostic plot of Logio([OIII]/H/3) vs 
Logio([NII]/Ha). AGN are defined to lie above the solid 
line, Hll-region like galaxies below the dashed line, and compos- 
ite galaxies between these boundaries. The filled circles indicate 
the emission line measurements for type 2 objects in the 2MASS 
sample, the stars represent narrow-line emission from the type 
1—1.9 objects, the triangles represent the objects in the 2MASS 
sample, which have fainter emission than the rest of the sample, 
and the small points are from SDSS-DR 8 data. 



rate narrow H/3 emission line fluxes in many of the objects, 
related to the de-blending of the emission of the BLR and 
the NLR (see discussion in ^ 3.3). 

Based on the diagnostic diagrams, all of the objects that 
could be plotted fall in the AGN part of the diagrams with 
the exceptions of J0411-01 and JlOOl+41, which fall in the 
composite/HII region part of the diagrams, suggesting a ma- 
jor contribution from stellar photoionization. We therefore 
exclude the latter two objects from any subsequent analysis 
involving AGN properties. 

In addition to emission line diagnostic ratios plotted in 
Figures [9] - the presence of FHILs in the spectra can in- 
dicate unambiguously the presence of an AGN. These lines 
include [NeV]AA3345,3426, [FeVII]A6086 and [FeX]A6374 
which, due to their high ionisation potentials (>97 eV), are 
not generally associated with LINER nuclei or starbursts. 
We find that 23 of the objects in our sample (85%) show 
evidence for at least one FHIL, suggesting significant AGN 
activity. 

Another indicator of the presence of an AGN is the 
detection of broad line emission from the BLR. Table [5] in- 
dicates whether broad emission lines are detected for each 
of the 2MASS objects. We assume that the BLR is detected 
if the Balmer recombination lines show significant (i.e. well 
above the noise) broad wings that cannot be accounted for 
by the [OIII] model, and have FWHM > 1000 km s'^ . 16 



Figure 10. Diagnostic plot of Logio([OIII]/H/3) vs 
Logio([OI]/HQ:). AGN are defined to lie above the solid 
line, with the dashed line separating Seyfert-type objects (above 
the line), and LINER-type objects (below). Hll-region like 
galaxies fall below the solid line. The symbols are the same as 
Figure |9] 



out of 27 objects in our sample (59%) show the presence of 
definite BLR emission. In addition, tentative evidence for 
weak Ha broad line emission in J1131+16 was presented in 
[Rose et al.| ( |20lT| . The broad line component in J1131-I-16 
is consistent with the presence of a scattered (Antonucci & 
'Miller 19851, or a directly observed, but weak, broad line 



region (BLR) component. The rest of the type 2 objects 
had no such emission. Figure [12] presents the velocity widths 
(FWHM) of the BLR components as measured from single 
Gaussian fits to the broad Balmer line^ plotted against 
hyoiii] for both the 2MASS sample and PG sample. Where 
possible, the broad and narrow-line emission were fitted sep- 
arately in order to obtain the BLR FWHM; such cases are 
indicated by the filled circles. Where the broad and narrow 
components could not be confidently separated we use the 
velocity widths of single Gaussian fits to the whole profiles; 
these cases are indicated by the black triangles. The me- 
dian BLR velocity width (FWHM) for the 2MASS sample 
is 2100±400 km s~^, which is lower than that of the PG sam- 
ple (4500±300 km s"^). However, the 2MASS median may 
be affected by the inclusion of NLR emission in the objects 
for which we could not confidently separate the emission 
from the BLR and the NLR. The median velocity width for 



^ Note that the FWHM of the Balmer lines are based on fits to 
the H,9 in all cases except of J1006-I-41 and J1040-I-59, which do 
not have any detectable broad emission in the H/3 emission line, 
but clear BLR Ho components. 



On the nature of the red, 2MASS selected AGN in the local Universe I: an optical spectroscopic study 




Figure 11. Diagnostic plot of Logio([OIII]/H/J) vs 
Logio([SII]/HQ:). AGN are defined to lie above the solid 
line, with the dashed line separating Seyfert-type objects (above 
the line), and LINER-type objects (below). Hll-region like 
galaxies fall below the solid line. The symbols are the same as 
Figure |9] 



the 2MASS objects in which the broad and narrow compo- 
nents could be confidently separated is 4400±1100 1cm s~^ 
— comparable with the PG quasars. 

All the results from this section are summarised in the 
truth table of Table [s] Overall, 25/27 (~93%) objects show 
some sign of AGN activity; and only 2 objects — J0411-01 
and JlOOl-l-41 — show no clear evidence of AGN activity. 
Given that it is not clear whether the latter two objects are 
powered by AGN activity, we do not include them in the 
analysis in the rest of the paper. In terms of a more detailed 
classification of the spectra of the 2MASS objects for which 
we detect significant AGN activity, 59% qualify as Seyfert 

1 galaxies, 30% as Seyfert 2 galaxies, and 4% have weak 
broad lines placing them between the Seyfert 1 and Seyfert 

2 classification. 



4-1.3 Narrow-line Seyfert 1 galaxies 

Narrow-line Seyfert 1 galaxies (NLSl) are a sub-class of 
Seyfert galaxies with BLR Balmer emission lines that are 
relatively narrow (FWHMh;^ < 2000 km s~^:|Osterbrock fc 



Pogge||1985| [Veron-Cetty, Veron fc Gongalves||2001[ ). There 
are several hypotheses surrounding the nature of these ob- 
jects. One prevailing idea is that NLSl are AGN with lower 
than typical black hole masses that are undergoing a phase 
of rapid accretion ( Pounds, Done fc Osbourne|1995 1. In this 
scenario the NLSl fall below, but are growing towards, the 



Figure 12. Velocity width (FWHM) of the BLR components 
plotted against Ij[oiii]- 2MASS objects for which broad 

emission can be confidently separated from the narrow line emis- 
sion are indicated by the filled circles, whereas the 2MASS ob- 
jects for which the broad and narrow emission could not be con- 
fidently separated are indicated by the filled triangles. The PG 
quasar population are indicated by the unfilled squares. All the 
line measurements are based on single Gaussian fits. 



they accrete more mass, they will appear as typical Seyfert 
objects ( |Mathur|200"0l ). 

If the red 2MASS objects are young, recently triggered 
AGN, then we might expect there to be a relatively large 
fraction of NLSl objects in the sample. We checked our sam- 
ple for NLSl based on single Gaussian fits to their H/3 lines, 
to find objects with H/3 or Ha emission line widths that are 
below the NLSl limit (FWHM < 2000 km s"^), but greater 
than the widths of the [OIII] lines. This is consistent with 
the fitting procedures of [Osterbrock fc Pogge| ( |1985^ and 
Goodrich (19891. We found that 4 objects with adequate 



single Gaussian fits to H/3 have velocity widths consistent 
with them being NLSl: J0248-f 14, J0422-18, J0435-06 and 
J1006-f-41. These objects have H/3 velocity widths (FWHM) 
of 1550±20 km s'S 1450±120 km s'^ , 1780±70 km s"^ and 
1170±70 km s"^ respectively. 

Figure [13] shows the single Gaussian fits to the H/3 emis- 
sion lines in these objects. Apart from the line widths, there 
are two other key properties that are found to be associated 
with NLSl: [OIII] emission lines that are weak relative to the 
broad H/3 lines ([OIII]/H/3tota!<3) and the presence of strong 
Fell emission (FeIl457o/H/3tota!>0.5: Osterbrock & Pogge 



1985 



|Goodrich| [19891 |Veron-Getty, Veron fc Gongalve" 



M-(j correlation (Ferrarese fc Merritt 20001; eventually, as 



2001|. We find that all 4 objects with FWHMjj,3<2000 km 
s~^ in our sample also posses these typical NLSl properties. 

Of the AGN in the 2MASS sample, only ~16% are NLSl 
objects, or ~24% as a fraction of all the broad emission-line 
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Q0248 - Aperture 1 




Q0422-1854 - Aperture 1 




Q0435-0635 - Aperture 1 




Table 5. Truth table presenting the evidence for AGN activ- 
ity in the 2MASS sample objects. A true result for BLR indi- 
cates the detection of a broad component for either/both Ha 
and H/3. A true result for the FHIL indicates the positive de- 
tection of a high ionisation species (>54.4 eV). The [Nil] col- 
umn shows the results from Figure [9] the [OI] column shows Fig- 
ure |10| and the [SII] column shows Figure T=true, F=false, 
P=tentative evidence for a BLR, A=AGN, S=Seyfert type, 
H=starburst galaxy, C=composite galaxy, D=unmeasured due to 
the degeneracy between the broad and narrow Balmer emission, 
and X=unmeasured because an appropriate emission line has not 
been detected. 



Name 


quasar 


BLR 


FHIL 


[Nil] 


[OI] 


[SII] 


0221-1-13 


F 


F 


T 


A 


S 


s 


0248 -M4 


F 


T 


T 


D 


D 


D 


0306-05 


F 


F 


T 


D 


D 


D 


0312+07 


F 


T 


T 


D 


D 


D 


0400-f05 


T 


T 


T 


A 


S 


S 


0409-f08 


F 


T 


T 


A 


S 


s 


0411-01 


F 


F 


F 


C 


X 


H 


0422-18 


F 


T 


T 


D 


D 


D 


0435-06 


F 


T 


F 


D 


D 


D 


0447-16 


F 


T 


T 


D 


D 


D 


0504-19 


F 


F 


T 


A 


S 


S 


0910-1-33 


F 


F 


T 


A 
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S 


1001+41 


F 


F 


F 


C 


H 


H 


1006+41 


F 


T 


T 


D 


D 


D 


1014+19 


F 


T 


T 


D 


D 


D 


1040+59 


F 


T 


F 


D 


D 


D 


1057-13 


F 


F 


T 


A 


S 


S 


1127+24 


F 


T 


T 


D 


D 


D 


1131+16 


T 


pa 


T 


A 


S 


S 


1158-30 


F 


F 


T 


A 


S 


S 


1212-14 


F 


T 


T 


A 


X 


s 


1321+13 


F 


F 


T 
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1323-02 


F 


T 
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s 
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1338-04 


T 


T 
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1407+42 


T 


F 
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1448+44 


F 


T 


T 


D 


D 


D 


2121-17 


T 


T 
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D 


D 


D 



See jRose et al.| l (2011| for full dcscripti( 



Q1006+4104 - Aperture 1 




Figure 13. Gaussian fits to the H/3 emission lines of the NLSl 
candidates in the 2MASS sample. Prom top to bottom: J0248+14, 
J0422-18, J0435-06 and J1006+41. The overall fits to the emis- 
sion lines are adequate and therefore we use the velocity widths 
(FWHM) to characterise these objects. 



objects in our 2MASS sample. This is similar to the fraction 
of NLSls found in other samples of broad-line AGN: ~24% 
for the unobscur ed type 1 AGN (|Jin et al.||2012a| l; 11% for 
the PG quasars ( |Hao et al.|[2005|); and ^5% o f aU QSOs 
in SDSS DR3 AGN catalogue JZhou et al.|2006[ ). Therefore 
we conclude that the 2MASS sample does not contain an 
unusually high proportion of NLSl objects. 



4.2 Balmer decrements and reddening 

One possible explanation for the red colours of the 2MASS 
objects at near-IR wavelengths is that their AGN suffer 
an unusual amount of dust extinction compared with other 
samples of AGN ( e.g. Glikman et al.| [2004). In order to in- 
vestigate this possibility, we have calculated Balmer recom- 
bination line decrements for the 2MASS, PG/SDSS and |Jin| 
|et al.||2"012a,b,c| samples to determine whether the 2MASS 
objects have higher levels of dust extinction than 'typical' 
blue quasars. 

Ideally, Balmer decrements (Fj/a/Fn^f) would be de- 
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termined for the narrow and broad line components of each 
object separately. However, some objects with a detectable 
BLR have relatively narrow broad lines (e.g. J1006+41), 
which makes it difficult to separate the components at the 
resolution of our spectra (see § 3.3). Therefore we have 
taken two approaches. For the first approach we used the 
measured total (broad and narrow) fluxes of the Balmer 
emission lines. The second approach was to separate the 
broad and narrow emission components of the blends, where 
the components could be confidently separated. We only 
used the total (broad + narrow) Balmer flux measurements 
for the Jin et al. (2012a I sample, because we were uncer- 



tain whether the broad and narrow emission components 
could be confidently separated in that sample. In addition, 
when we consider Balmer decrements for the total flux of the 
Balmer lines, we only use the data for the objects with a de- 
tectable BLR, and do not consider the Balmer decrements 
for the type 2 objects. Balmer decrements for the 2MASS 
objects are given in Table |6] 

Figure [14] shows histograms of the Balmer decrements 
{P Ha/'P Hfj) for the combined broad and narrow Balmer 
emission of the 2MASS sample (broad line objects only), 
the PG/SDSS sample and the unobscured type 1 AGN. 

It is clear from this figure that, on average, the 2MASS 
sample objects have higher PHa/PHfi ratios than the PG 
quasars and unobscured type 1 AGN, implying that there 
is more dust extinction in the 2MASS sample. There is 
also a significant fraction of objects (5/19) with Balmer 
decrements that are higher than the most reddened PG 
quasars. This is reflected in the median "PHa/PHfi ratios: 
F//c«/F//;3=4.9±0.5, 3.2±0.4, and 2.8±0.1 for the 2MASS, 
PG/SDSS and unobscured type 1 objects respectively. 

To put this in context, the Balmer decrements for the 
2MASS sample suggest reddening in the range 0.0 ^ E(B-V) 
^ 1.2, with a median E(B-V)=0.52±0.05, assuming a Case 
B Ha/H/3 ratio of 3.1 ( [GaskeU fc FerlandlpSif . However, 
although the 2MASS sample displays a significant degree of 
dust reddening, there are several objects (~25% of all the 
broad line objects) which show very little (if any) significant 
BLR reddening. These objects include J0312-I-07, J0409-I-08, 
J0435-06, J1338-04 and J2121-17. This is interesting, be- 
cause it highlights the diverse nature of the 2MASS sample. 
We further note that three of the objects with low Balmer 
decrements — J0312-I-07, J0435-06 and J1338-04 — have J- 
Ks colours just below the J-Ks=2.0 limit. Also, J1040-)-59, 
an object with one of the highest Balmer decrements, has 
one of the reddest J-Ks colours in the sample (J-Ks=3.02). 
However, overall we do not find a statistically significant 
correlation between J-Kg colours and Balmer decrements 
for the 2MASS sample. 

We have performed K-S tests ( Peacock|1983 Fasano & 



Franceschini 19871 to check the significance of any differ- 



ences in the distribution of the Balmer decrements between 
the samples. We find that we can reject the null hypothe- 
sis that the 2MASS sample is drawn from the same parent 
population as the PG quasars and unobscured type 1 AGN 
at the 2a and >3(t levels of significance respectively (the 
corresponding p values are 0.007 and <0.001). 

In the cases where it was possible to clearly separate the 
BLR and NLR components, we find that the Balmer decre- 
ments obtained for the BLR components alone are not sig- 
nificantly different from those derived from the total (BLR 




10 12 



Balmer Decrement 




4 6 

Balmer Decrement 



10 12 



Unobscured type 1 AGN 



Balmer Decrement 

Figure 14. Comparison of the total flux Balmer decrements 
{P Ha/P Hp)- The 2MASS population measurements (broad line 
objects only) are indicated in black (top), the PG quasar popula- 
tion are shown in grey (middle) and the unobscured type 1 AGN 
are given in white (bottom). Note that the Balmer decrements 
are only presented for those objects with broad Balmer emission. 
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Table 6. The Ha/H/3 Balmer decrements measured for the 
2MASS sample. The second column presents the Balmer decre- 
ments for the total (BLR + NLR) fluxes. The third and fourth 
columns present the NLR and BLR Balmer decrements respec- 
tively. 
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J0447-16 


8.74 it 1.48 


_ 


_ 


J0504-19 




3.83 it 0.82 




J0910-I-33 




2.15 it 0.62 




J1006-I-41'' 


5.96 it 0.28 




>5.00 


J1014-I-19 


5.57 it 0.79 






J1040-I-59'' 


12.14 it 0.57 




>6.67 


J1057-13 




8.05 it 0.67 




J1127-I-24 


4.89 it 0.46 


3.25 it 0.46 


5.38 it 0.60 


J1131-l-16'= 




5.00 it 0.17 




J1158-30 




10.98 it 0.48 




J1212-14 


4.79 it 0.93 


5.07 it 0.93 


4.71 it 0.75 


J1321-I-13 




5.67 it 0.52 




J1323-02 


4.75 it 0.87 






J1338-04 


3.75 it 0.69 






J1407-I-42 




6.79 it 0.68 




J1448-f44 


4.90 it 1.21 






J2124-17 


2.15 it 1.19 








Figure 15. The Balmer decrement (Fhol/Php) versus rcdshift, 
z, for the NLRs of the 2MASS sample. The 2MASS population 
measurements are indicated by the filled circles, the PG/SDSS 
quasars are indicated by unfilled squares and the unobscured type 
1 objects are indicated by asterisks. 



" Sky line features cut through the region of the H/S emission 
line. Results are based on an upper limit measurements for 
the H/9 emission line, but are presented as a lower limit on 
the Balmer decrement. 

BLR Balmer decrement is based on an upper limit of the H/3 
flux and therefore is a lower limit of the Balmer decrement. 

The emission flux from Ha is boosted by coUisional excita- 
tion in this object jRose et al.|20lTl l. 

+ NLR) fluxes (see Table [6]|. Therefore we do not consider 
the pure BLR Balmer decrements further in this paper. 

When considering the range of reddening found for our 
2MASS AGN, our results are consistent with the idea that 
dust reddening plays a significant role in causing the red 
near-IR colours of the 2MASS objects (see ^ 6.1). 



4-2.1 Reddening in the NLR 

It is also interesting to consider the NLR Balmer decrements 
for the red 2MASS AGN in which the NLR and BLR compo- 
nents could be clearly separated. The results are presented 
in Figure [TS] We find that the NLR Balmer decrements cover 
a similar range to those measured for the BLR (and BLR 
+ NLR combined), with some objects showing evidence for 
a high degree of reddening, but others showing no signifi- 
cant reddening. We have compared the distributions of the 
narrow emission line Balmer decrements for the 2MASS ob- 
jects and unobscured type 1 AGN using a K-S test and find 
that we cannot reject the null hypothesis that these samples 



are drawn from the same parent population, because the test 
gives a high p value: 0.099. Note that small number statistics 
preclude a detailed comparison of the NLR Balmer decre- 
ments between the 2MASS and PG/SDSS samples. Also, in 
the case of the NLR it may not be valid to assume Case 
B recombination for the intrinsic reddened Balmer decre- 
ment: our detailed analysis of J113H-16 ( Rose et al.||2011 l 
demonstrates that coUisional excitation of the Ha emission 
line in the partially ionised zone can significantly enhance 
the Ha/H/3 ratio if the gas densities are high (e.g. if part of 
the NLR emission originates in the circum- nuclear torus). 



4.3 Outflows 

If the red 2MASS objects truly represent young AGN en- 
shrouded in a cocoon of dust and gas, then we might expect 
them to show strong evidence for outfiows in their spec- 
tra, when compared to other types of active galaxies. For 
example, in models of major galaxy mergers the AGN are 
triggered as gas and dust is driven by the tidal torques asso- 
ciated with the mergers into the nuclear regions of the galax- 
ies; the AGN accrete at high Eddington rates, and then drive 
powerful outfiows, blowing away the enshrouding dust and 
gas (|lfo^ ms et al.|2005||Naray anan et al.|2006[ [Hopkins et| 
al.|2006[ ). 

Figure [16] shows a histogram of the rest-frame velocity 
shifts of the broad [OIII] components (where detected) rela- 
tive to those of the narrower components, as calculated from 
the parameters given in Table [2] The 2MASS objects are 
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Figure 16. Histogram of the shift of the broad [OIII]A5007 com- 
ponents relative to the narrow components, as measured using 
DIPSO. The 2MASS measurements are indicated by the bla<;k bars 

(top), the PG quasar poijulation by the grey bars (middle) and 
the unobscured type 1 AGN by the open bars (bottom). 
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Figure 17. Histogram of the instrumentally corrected velocity 
widths (FWHM, km s'^) of the broad [Oni]A5007 component, as 
measured using dipso. The 2MASS measurements are indicated 

by the black bars (top), the PG quasar population by the grey 
bars (middle) and the unobscured type 1 AGN by the open bars 
(bottom). 
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Figure 18. Outflow velocity (km s ^) plotted against the ve- 
locity width at FWHM (km s^i) of the broad [OIII]A5007 com- 
ponents as measured by DIPSO. The 2MASS population measure- 
ments are indicated by the filled circles, the PG quasar population 
are indicated by the unfilled squares and the unobscured type 1 
AGN are indicated by asterisks. 



nents. As with the velocity shifts, the widths (FWHM) of 
all three samples have a comparable range. Again, we have 
tested for any significant differences between the distribution 
of velocity widths (FWffM) of the samples using a f D K-S 
test. Between the 2MASS objects and PG quasars we find 
a p value=0.385, and a p value=0.089 for the unobscured 
type 1 AGN. Therefore we cannot reject the null hypothesis 
that the samples are drawn from the same parent popula- 
tion when comparing the linewidths of the PG quasars and 
unobscured type 1 AGN with the 2MASS sample. 

Figure [is] shows the outflow velocity plotted against the 
velocity width (FWHM) of the broad [OIII]A5007 compo- 
nent. We present this Figure to further illustrate that there 
are no significant differences between the [OHI] emission line 
kinematics of the samples. We have also tested for any signif- 
icant difference between the populations using a 2D K-S test 
( Peacock|1983 Fasano fc Franceschini|f987) . Again, we find 
no evidence for significant differences between the samples 
(p > 0.05). 

One important caveat when considering these results is 
that the emission line kinematics may be affected by the ori- 
entation of the AGN. One possibility is that our line of sight 
to the type 2 AGN is at a larger angle to the torus axis than 
type f AGN (as suggested by the orientation-based unified 
schemes). Therefore due to projection effects the measured 
blueshifts of the broad components might be lower in the 
type 2 than the type 1 objects. Alternatively, it is possible 
that the torus obscures the narrow line outfiows from our di- 
rect line of sight in the type 2 objects. However, comparisons 
between the NLR kinematics of type 1 and 2 AGN show no 
clear differences (e.g. |Whittle|1985[ [Nelson fc Whittle|1995 1. 



plotted along with the PG quasars and the unobscured type 
1 AGN for comparison. For the unobscured type f AGN, 



we used the online data tables from Jin et al. (20f2al to 



determine the outflow properties. For the 2MASS objects, 
PG quasars and unobscured type f AGN, we assume that 
the narrow component of the [OHI] emission represents the 
host galaxy rest frame emission, and the broader-line emis- 



sion component the outflow (Holt, Tadhunter & Morganti 



2008) 



From Figure [16] it is clear that there are no signiflcant 
differences between the distributions of the velocity shifts of 
the samples and, the broad components are predominantly 
blue-shifted in all three samples, with only a few exceptions. 
We flnd the median velocity shifts of -240±80 km s~^, 
-230±50 km s'^ , and -160±30 km s"^ for the 2MASS, PG 
and unobscured type I samples respectively. We have tested 
the signiflcance of any difference between the distribution of 
velocity shifts of the samples using a ID K-S test. Between 
the 2MASS objects and PG quasars we find a p value=0.731, 
and between the 2MASS objects and unobscured type 1 ob- 
jects we find p value=0.230. Therefore we cannot confidently 
reject the null hypothesis that these samples are drawn from 
the same parent population. Note that, although there are 
no significant differences between the three samples, the 
most extreme object in all three samples in terms of its shift 
is a red 2MASS object: J1158-30. This object might be a 
good candidate for a young AGN. 

In addition, in Figure[l7]we plot a histogram of the rest- 
frame velocity widths (FWHM) of the broad [OHI] compo- 



5 BLACK HOLE AND ACCRETION 
PROPERTIES 

Further insights into the nature of the red 2MASS objects 
can be obtained by studying the SMBH and accretion prop- 
erties of the sample. Young, dust enshrouded objects are 
expected to have higher than average accretion rates, be- 
cause there is more gas available for accretion (e.g. see pl| 
Matte o, Sprin gel fc Hern quist "2005[) . 

To determine the accretion properties, we must first de- 
termine the black hole masses of the objects. In this section 
we use the broad Balmer emission line properties to deter- 
mine virial black hole masses and Eddington ratios for the 
2MASS sample, the PG sample, and the unobscured type 
1 AGI^Q Any significant differences in the Eddington ra- 
tios between the samples could indicate differences in their 
evolutionary stage. 

To determine virial black hole masses, we have used the 
relationships between the Balmer emission and black hole 



^ For consistency, we use the H/3 luminosities and velocity widths 
(FWHM) provided in |Jin et al.|2012a| to determine the black hole 
masses, bolometric and Eddington luminosities of the unobscured 
type 1 AGN sample (L^^ and FWHMj:/^ are flux weighted means 
of the intermediate and broad components of the H/3 emission 
lines). Note that, the data used from [jin et al.|2012a| is not cor- 
rected for extinction. This is reasonable as the sample is selected 
to be unreddened. 
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mass from Greene & Ho ( 2005 1 , which are based on corre- 



lations between the AGN bolometric luminosities {Lbol), 
Balmer lino luminosities, and line widths (FWHM) (see 
Greene & Ho 2005, equations 6 and 7). We use equation 



7 for most cases where we can use the broad H/3 emis- 
sion lines, minimising any degeneracies due to blending with 
other emission lines (e.g. [Nil] with Ha). In cases where the 
H/3 broad emission is not detectable, we used the Ha broad 
emission with equation 6 instead. 

We fitted the broad H/3 emission with the minimum 
number of Gaussian components required to obtain an ad- 
equate fit to the emission from the BLR, and took a fiux 
weighted mean of the FWHM, consistent with the method 



of Greene & Ho (20051. Where possible, we corrected the 



BLR fluxes for extinction, and in cases where we could not 
confidently separate the broad and narrow emission, we used 
the total Balmer decrement measurements for the extinc- 
tion correction. Note that, while our approach is consistent 



with that outlined in Greene & Ho ( 2005 1 when measuring 



the broad H/3 emission, problems with the separation of the 
broad and narrow emission in the objects where we could 
not separate the components may potentially have an affect 
on the line width measurements. Also, for the cases where 
we estimated an upper limit on the broad H/? flux, we used 
the lower limiting estimate of the extinction to flnd a lower 
limit on the SMBH mass, and therefore an upper limit on 
the Eddington ratio. The 2MASS objects have virial SMBH 
masses in the range 7.2< logio[Ms_ff]<9.2 Mq. 

As well as making virial SMBH estimates based on 
the broad emission line properties, for comparison we have 
also used the absolute R-band magnitudes (M^) of the host 
galaxies of the 2MASS sample taken from |Hutchings et al.| 
(20031. These estimates have the advantage of being free 



from any degeneracies caused by the blending of the emis- 
sion lines from the BLR and the NLRs of the sample. In 
addition, masses can be estimated for all the sources, rather 
than just the sources with an observable BLR. To make these 
estimates we used the relation given in |McLure fc Dunlop| 
(20011, which takes advantage of the correlation between 



black hole mass and host galaxy bulge luminosity, to deter- 
mine the SMBH masses. 

Figure [l9] shows the SMBH mass estimates from the 
R-band magnitudes plotted against the masses calculated 
using the properties of the BLR. In addition, we have plot- 
ted a line showing a one-to-one relationship between the 
black hole masses. While some objects show good agree- 
ment between the two estimates, there is a large scatter 
which is reflected in the uncertainty of the median ratio 
{M^fria? /^R-ba"d)- 0.81±0.81. lu cascs whcrc the virial 
mass estimates are lower than the R-band estimates, orien- 
tation/geometric and/or reddening effects could potentially 
lead to the Balmer fluxes, velocity widths and SMBH masses 
being underestimated. Alternatively the R-band measure- 
ments may be contaminated by light from the disks of the 
host galaxies, leading to over-estimates of the SMBH masses. 
Certainly, at least one object in our sample - J1131-I-16 - 



has a strong disk component (Rose et al. 20111. However, 



in cases where the virial masses are higher, the flux of the 
BLR may be overestimated because of an inaccurate mea- 
surement in the broad Ha flux due to blending with the 
narrow Ha-|-[NII] lines, which could lead to an overestima- 
tion of the level of extinction. 
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Figure 19. A comparison of the SMBH mass estimates for the 
2MASS AGN. We plot the mass estimates based on the the R- 
band magnitudes against the virial mass estimates based on the 
broad Balmer emission lines. In addition, we have plotted a line 
showing a one to one relationship for the black hole masses. The 
mass estimates obtained from the R-band magnitudes are higher 
on average than those obtained from the BLR measurements. 



Figure [20] shows the virial black hole masses plotted 
against Ijyoiii] for the 2MASS and comparison samples, 
where the uncertainties are only plotted for the 2MASS ob- 
jects for clarity. We do not plot the R-band black hole mass 
estimates because we do not have R-band magnitudes for 
the comparison samples, and we want to present the compar- 
isons between mass estimates which have been determined 
via the same method. In addition. Table [7] shows the indi- 
vidual SMBH masses of the 2MASS sample. Overall, from 
this figure it is clear that the black hole masses of the PG 
quasars (unfilled squares) and unobscured type 1 AGN (as- 
terisks) are generally lower than the 2MASS objects (filled 
symbols) of equivalent [OIII] luminosity. Potentially, this dif- 
ference can be explained in terms of extinction effects: in the 
cases of the few objects with virial mass estimates for which 
we have been able to correct hyoiii] for extinction in the 
NLR, the extinction correction moves the points onto the 
main correlation deflned by the PG quasars and unobscured 
type 1 objects. 

Alternatively, the black hole masses of the 2MASS ob- 
jects could be genuinely higher than those of the PG quasars 
and unobscured type 1 AGN of similar intrinsic power. In 
this context we note that the recent study of |Canalizo et aX\ 
\2Q12\ measured black hole masses and host galaxy proper- 
ties for 29 red 2MASS AGN with 0.17 < z < 0.37, flnding 
that the majority of red 2MASS AGN also have significantly 
more massive black hole masses than the masses of black 
holes in local AGN. However, we remain cautious in our in- 
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Table 7. Virial SMBH masses and Eddington ratio results for the 2MASS sample. Values presented in brackets are for those objects 
where hgoi^ has been corrected for dust extinction. 



Name Virial (^bol I'^edd) 

Logio[MsMBff] (M0) Logio[MsA/SH] (M0) 



J0221+13 


8.19 


± 


0.43 




- 




- 


J0248+14 


8.07 


± 


0.44 


7.46 


± 


0.09 


0.020 ± 0.005 


J0306-05 


8.25 


± 


0.43 




- 




- 


J0312+07 


8.07 


± 


0.44 


8.29 


± 


0.32 


0.009 ± 0.001 


J0400+05'' 


7.95 


± 


0.44 


9.14 


± 


0.08 


0.056(0.585) ± 0.009(0.125) 


J0409+07'' 


7.58 


± 


0.46 


7.70 


± 


0.05 


0.204(0.901) ± 0.025(0.112) 


J0422-18 


8.01 


± 


0.44 


7.92 


± 


0.11 


0.079 ± 0.024 


J0435-06 


8.68 


± 


0.41 


7.87 


± 


0.26 


0.065 ± 0.003 


J0447-16 


8.19 


± 


0.43 


8.10 


± 


0.11 


0.023 ± 0.007 


J0504-19 


8.01 


± 


0.44 










J0910+33 


7.83 




0.45 










J1006+41'' 


7.52 




0.46 


7.26 




0.10 


0.302 ± 0.105 


J1014+19 


8.44 


± 


0.42 


8.80 


± 


0.13 


0.041 ± 0.014 


J1040+59'' 


7.52 


± 


0.46 


8.40 


± 


0.11 


0.021 ± 0.007 


J1057-13 


7.77 


± 


0.45 










J1127+24° 


7.64 


± 


0.45 


7.80 




0.10 


0.036(0.041) ± 0.019(0.0010) 


J1131+16 


8.31 


± 


0.43 










J1158-30 


8.25 


± 


0.43 










J1212-14'^ 


7.77 


± 


0.45 


7.16 




0.18 


0.037(0.172) ± 0.020(0.090) 


J1321+13 


7.52 


± 


0.46 










J1323-02 


8.38 


± 


0.42 


8.58 


± 


0.16 


0.057 ± 0.026 


J1338-04 


8.50 


± 


0.42 


7.96 


± 


0.26 


0.383 ± 0.025 


J1407+42 


7.95 


± 


0.44 










J1448+44 


8.44 


± 


0.42 


7.75 




0.21 


0.222 ± 0.137 


J2124-17 


8.68 


± 


0.41 


8.27 


± 


0.04 


0.802 ± 0.068 



" The NLR of this object has been corrected for dust extinction. 

^ ^SMBH E'H'l therefore a^re calculated using broad Ho emission as opposed to H/3. 



terpretation of the derived black hole masses for our 2MASS 
sample because of the uncertainty surrounding the extinc- 
tion corrections. 

In terms of estimating the Eddington ratio 
(L_bol/L_b_dd), we used hyoiii] to calculate the bolo- 
metric luminosities of the AGN, rather than the AGN 
continuum, because the [OIII] emission line is well detected 
in all the objects, whereas the AGN continuum emission 
is subject to reddening effects and potential contamination 
by the host stellar continuum. Lbol was determined from 



L[o/jj] usin g the relationship Lsol=3500L[o//7] (Heckman 
et al.|2004| |. 



Figure [21] shows the Eddington ratios plotted against 
LsOL- Once again the uncertainties are only plotted for the 
2MASS objects for clarity. We present Eddington ratios for 
all possible objects in the 2MASS sample without correcting 
[OIII] for dust extinction (filled circles) and, where possible, 
we also show the Eddington ratios following the correction 
of the [OIII] luminosity for dust extinction (filled triangles), 
connecting both sets of estimates with a straight line. It is 
clear that, for the 4 objects in which we can measure red- 
dening accurately, reddening has a significant effect on the 
Eddington ratio. There is an obvious, large scatter in Ed- 
dington ratios for all the samples: the 2MASS objects have 
ratios in the range: 0.009 < Lbol /I^edd < 0.91. The me- 
dian ratio for the 2MASS sample is 0.05±0.05, without cor- 
recting [OIII] for reddening and 0.08±0.04 in the redden- 
ing corrected case. In comparison, the PG quasar sample 



(unfilled squares) and the unobscured type 1 AGN sample 
(asterisks) have median ratios of 0.22±0.12 and 0.23±0.14 
respectively; the large uncertainties reflect the large scatter 
in Eddington ratio. These ratios are all consistent within 
2a of the uncertainties. Note that, if we could correct for 
NLR reddening in all the objects, the median ratio for the 
2MASS objects would likely increase to a value closer to 
that measured for the PG quasar and unobscured type 1 
AGN samples. Clearly, there is no evidence that the 2MASS 
objects have Eddington ratios that are significantly higher 
than those of the comparison samples. 



6 DISCUSSION 

6.1 Moderately Reddened Objects 

A key result from this study is that many of the red 2MASS 
AGN have moderate amounts of dust extinction when com- 
pared to the comparison samples (0.2 < E(B-V) < 1.2); 
the reddening values we measure are consistent with those 
determined for other samples of red 2MASS AGN (0.1 < 
E(B-V) < 3.2: Glikman et al. 2004, 2007, 2011 & 2012; Shi 
et al. 2007; Urrutia et al. 2008, 2009 & 2012; Georgakakis et 
al. 2009; Kuraszkiewicz et al. 2009a,b; Canalizo et al. 2012). 
Therefore it is possible that the 2MASS AGN have red J-Ks 
colours because dust extinction reduces the J-band flux rel- 
ative to the K-band flux such that the Fk/Fj flux ratio is 
higher than a typical AGN. In Paper II we show that for an 
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Figure 20. Virial black hole mass (Mq) plotted against l^yoiii] 
( W) . The filled circles and triangles represent the 2M ASS sample. 
Circles represent black hole masses estimated using H/3 emission 

and filled squares arc based on the Ho emission. Where possible, 
we correct ^[oiii] for reddening, these estimates are represented 
by filled triangles and are joined to their uncorrected counterparts 
by a straight line. The unfilled squares represent the black hole 
mass estimates of the PG sample and the asterisks indicate the 
type 1 AGN estimates. 



Figure 21. Eddington ratio plotted against Ijbol (erg s~^) 
based on measurements of L[o///] E^^d virial black hole mass esti- 
mates. Filled circles represent the 2MASS objects where l^yoiii] 
has not been corrected for dust extinction. For objects which 
could be corrected for dust extinction, we represent them with 
filled triangles and join them to their respective estimates where 
no correction has been made. The unfilled squares represent the 
PG quasars and the asterisks represent the unobscured type 1 
objects. 



AGN with intrinsic near-IR colours that are typical of the 
PG quasars and the unobscured type 1 AGN (J-Ks~1.70), 
we require only a moderate level of reddening of E(B-V) > 
0.8 in order to give J-Ks > 2.0 — this level of extinction is 
measured for 29% of the 2MASS AGN with broad Balmer 
emission lines. 

Rather than the red near-IR colours being a result of 
dust extinction from a natal cocoon of dust, it is possible 
that the obscuring dust is located in the outer layers of the 
torus, and that we are observing these objects at an angle 
where our line of sight grazes the edge of the torus, leading 
to the observed red colours. However, such modest levels of 
reddening could also be produced by dust on a kpc-scale in 
the disks of the host galaxies. . 

Alternatively, rather than moderately extinguished 
AGN, it is possible that some of these objects are genuine 
type 2 objects in which the torus completely extinguishes 
the broad-linc AGN in the optical and shorter wavelength 
near-IR bands. In this case, the level of extinction is ex- 
pected to be much higher (A\/ > 10) than that indicated 
by the Balmer decrements in the NLR because, while the 
SMBH and BLR are completely obscured at optical wave- 
lengths by the torus, the NLR extends beyond the scale of 
the torus and is therefore relatively unaffected. In this case, 
due to the much higher levels of extinction caused by the 
torus, their J-Kg colours should be much redder than is ob- 



served. Indeed, from our work in Paper II, we find that A„ > 
5 would produce J-Ks > 3.00, a level only found for the type 
2 2MASS object J1307-I-23 (excluded from the main study 
of this paper), but not for any of the other type 2 objects 
in our sample. However, the integrated near-IR colours of 
the heavily extinguished type 2 objects could be made bluer 
if the host galaxies make a significant contribution in the 
near-IR, because the colours of typical unreddened stellar 
populations are relatively blue at such wavelengths. 

Finally, we note that dust extinction is not the only ex- 
planation for the red near-IR colours of the 2MASS objects. 
Some objects in our sample show relatively modest levels of 
extinction based on their BLR Balmer decrements, and have 
optical continuum shapes that appear similar to UV/optical 
selected AGN (e.g. J2124-17). Rather than moderate levels 
of dust extinction in such cases, a relatively large covering 
factor for the hot dust (T ~ 1500K) could also result in red 
J-Ks colours. Emission from hot dust in the torus becomes 
significant at wavelengths which coincide with the K-band. 
Therefore, a larger than average covering factor for this hot 
dust could increase the fiux measured in the K-band rela- 
tive to that of the J-band (likely dominated by accretion disk 
emission), thus producing the red J-Kg colours measured for 
the 2MASS sample. 

All of these possibilities will be explored further in Pa- 
per II, where we will model the near- to mid-IR SEDs of the 



26 M. Rose, C. N. Tadhunter, J. Holt, J. Rodriguez Zaurin 



2MASS objects with varying the level of extinction, dust 
covering factors and host galaxy contributions. 



6.2 Young, dust enshrouded objects 

There have been several studies of red 2MASS AGN in the 
past that have focussed on different wavelength regions of 
the SEDs of the population (e.g. X-ray, UV, optical or IR). 
In this section we outline the basic properties and findings of 
these studies, and compare them with the results obtained in 
this paper, in order to determine the nature of our red, low- 
z, 2MASS AGN sample. The basic properties of the samples 
used for the previous studies are summarised in Table [8] 

Many of these past studies conclude that the red 
2MASS quasars are young, dust-enshrouded, transitional 
objects. This conclusion is based on several lines of evi- 
dence, including a high rate of morphological disturbance 
of the host galaxies ( Urrutia, Lacy fc Becker]|2008 1, a rela- 
tively high rate of occurrence low ionization broad absorp- 



tion lines (LoBALs; Urrutia et al. 20091, evidence for en- 



hanced star formation rates compared to typical UV/optical 
selected AGN ( |Shi et al.|2007[|Georgakakis et al.|2009[ ), and 
high Eddington ratios for the gas being accreted by the cen- 
tral super-massive black holes l Urrutia et al.|2012 |. 

In contrast, although the results from our local 2MASS 
AGN sample agree with the these previous studies in the 
sense that we determine the same moderate levels of red- 
dening, we find no strong evidence that these are young 
quasars based on their [QUI] outflow properties, the inci- 
dence of NLSl, and the distribution of Eddington ratios. 
Indeed, the properties our 2MASS sample appear similar to 
those of the samples of UV/optical selected quasars in the 
local Universe. 

The apparent discrepancy between the conclusions of 
our study and the previous studies that identify the red 
2MASS objects as young, transition objects may be due to 
different selection criteria. In particular, the Glikman et al. 
(2004, 2007, 2011, 2012), Georgakakis (2009) and Urrutia et 
al. (2008, 2009, 2012) samples extend up to much higher red- 
shifts than our local sample, include optical colour selection 
criteria as well as near-IR colour selection, and some require 
the objects to have broad lines; in addition, Glikman et al. 
(2004, 2007, 2011, 2012) require detection of the objects in 
the FIRST radio survey. Such selection criteria favour more 
extreme objects that are more likely to be genuinely young 
and dust enshrouded. 

In this context it is notable that other studies of 2MASS 
AGN which are based on lower redshift samples and/or sam- 
ples selected without additional optical colour criteria (the 
bottom four entries in Table 8), agree with with the con- 
clusions of our investigation, albeit using different analysis 
techniques: 

- from X-ray to FIR SED studies, and PCA analysis 
of a population of red 2MASS AGN, |Kuraszkiewicz et al.| 
1 2009a,b|) conclude that the red 2MASS AGN are not ob- 



served at a younger evolutionary phase than other local sam- 
ples; 

- using I-band HST images to study the host galaxies of 
red 2MASS AGN, and comparing their findings to the PG 



quasars and a sample of ULIRGs, Marble et al. ( 2003 1 find 



no indication that red 2MASS AGN represent young, dust 
enshrouded, transitional objects; 

- by comparing optical broad-band polarimetry results of 
red 2M ASS AGN, PG quasars and BAL AGN, [Smith et al.| 
( 2002[ ) find that overall 2MASS AGN show relatively high 
levels of polarization (P > 3%), consistent with the idea 
that orientation/obscuration effects are responsible for the 
red NIR colours. 

These comparisons between the previous studies of the 
red 2MASS AGN emphasise the sensitivity of the conclu- 
sions to the precise selection criteria employed for the partic- 
ular samples. However, it is seems unlikely that the majority 
of low redshift 2MASS AGN that are selected solely on the 
basis of near-IR colour [[J ~ K) 2] are young, transitional 
objects. 



T CONCLUSION AND FUTURE WORK 

In this work we present a study of the optical spectra of 
a representative sample of 29 nearby (z < 0.28) 2MASS- 
selected AGN with red near-IR colours (J-Ks > 2.0). We 
compare the results from this sample with those obtained 
for comparison samples of PG quasars and unobscured type 
1 AGN. We find the following. 

• Overall the local 2MASS AGN population show sig- 
nificantly higher reddening than samples of UV/optical se- 
lected AGN. However, while some objects are highly red- 
dened (E(B-V) = 1.2), others show relatively little (if any) 
reddening. 

• The velocity shifts and widths of the [QUI] lines of the 
red 2MASS AGN are not significantly larger than those of 
the comparison samples. 

• The Eddington ratios of the 2MASS AGN and the in- 
cidence of NLSl type objects are comparable with those of 
'typical' AGN. 

These results do not support the view that the local 
red 2MASS AGN represent young quasars in a transitional 
phase, where the AGN are emerging from their natal cocoons 
of gas and dust. However, they are consistent with the idea 
that these objects are reddened either by obscuring dust in 
the outer layers of the circum-nuclear torii, or by dust in the 
kpc-scale disks of the host galaxies. 

In Paper II we will focus on the near-IR properties of 
this sample of red 2MASS AGN, and investigate the origin 
of their red J-Ks colours. 
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Table 8. Basic properties of red AGN samples. 



Sample 



Size 



Data used in Study 



Key Selection Criteria 



This paper 



Glikman et al. 2004, 2007, 
2011 & 2012 



Georgakakis et al. 2009 



Urrutia et al. 2008, 2009 
& 2012 

Shi et al. 2007 



Kuraszkiewicz 
et al. 2009a, 

Canalizo et al. 2012 



Marble et al. 2003 
Smith et al. 2002 



29 0.06 < z < 0.28 



Optical Spectra 



J -Kfl > 2.0, 



120 0.1 < z < 2.5 Optical & IR Spectra 



10 0.29 < z < 3.1 UV-FIR SED fitting 



13 0.1 < z < 1.0 



25 0.08 < z < 0.37 



Optical Spectra 



IR Spectra 



44 0.07 < z < 0.37 X-ray-FIR SED fitting, 

PCA analysis 

9 0.13 < z < 0.37 Optical/NIR spectroscopy, 

HST images 



29 0.13 < z < 0.6 
20 0.06 < z < 0.6 



HST images 

Optical broadband 
polarimetry 



selected from 



Hutchings et al. 



I 2003 I 



J-Kg > 1.7 & R-K ^ 4.0, 
selected by cross-correlating the 2MASS 
PSC with the FIRST radio catalogue 

J-Ks > 1.5 & R-K ^ 5.0, 
selected by cross-correlating the 2MASS 
PSC with SDSS DR3 spectra 

Selected from the 



Glikman et al. 



( 2007 1 sample 



J- K.^ > 2.0, 



selected from 



Cutri et al. 



1 2001 1 



J-Ks > 2.0 and B-Kg > 4.3 

J-Ks > 2.0 and < -25 

J-Ks > 2.0 and Uk < -23 
J-K5 > 2.0 and Mk < -23 
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the University of Massachusetts and the Infrared Process- 
ing and Analysis Center/California Institute of Technology, 
funded by the National Aeronautics and Space Adminis- 
tration and the National Science Foundation. This research 
has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the 
National Aeronautics and Space Administration. Funding 
for SDSS-III has been provided by the Alfred P. Sloan Foun- 
dation, the Participating Institutions, the National Science 
Foundation, and the U.S. Department of Energy Office of 
Science. 



REFERENCES 

Antonucci R.R.J. , Miller, J.S., 1985, ApJ, 297, 621 
Antonucci R., 1993, ARA&A, 31, 473 
Baldwin J., Philhps M., Terlevich R. 1981, PASP, 93, 5 
BeUamy M. J., Tadhunter C. N., Morganti R., Wills K. A., 

Hoh J., Taylor M. D., Watson C. A., 2003, MNRAS, 344, 

80 

Bian W., Gu Q., Zhao Y., Chao L., Cui Q., 2006, MNRAS, 
372, 876 

Boroson T. A., Green R. F., 1992, ApJS, 80, 109 
Bruzual G., Chariot S., 2003, MNRAS, 344, 1000 
Calzetti D., Armus L., Bohlin R. C, Kinney A. L., Koorn- 
neef J., Storchi-Bergmann Th., 2000, ApJ, 533, 682 



Canalizo G., Wold M., Hiner K. D., Lazarova M., Lacy M., 

Aylor K., 2012, ApJ, 760, 38 
Cutri R. M., Nelson B. O., Kirkpatrick J.D., Huchra J. P., 

Smith, P. S., 2001, ASPC, 232, 78 
Cutri R. M., Nelson B. O., Francis P. J., Smith, P. S., 2002, 

ASPC, 284, 127 
Dain S., Lousto C. O., Zlochower Y., 2008, PhRvD, 78, 

4039 

Dicken D., Tadhunter, C; Axon, D.; Morganti, R.; Inskip, 
K. J.; Hoh, J.; Gonzalez Delgado, R.; Groves, B., 2009, 
ApJ, 694, 268 

Dickson, R., Tadhunter, C, Shaw, M., Clark, N., Morganti, 

R., 1995, MNRAS, 273, 29 
Di Matteo T., Springel V., Hernquist, L., 2005, Nature, 

433, 604 

Fasano G., Franceschini A., 1987, MNRAS, 225, 155 
Ferrarese, F., Merritt, D., 2000, ApJ, 539, 9 
Francis P.J., Whiting M.T., Webster R.L., 1999, PASA, 17, 
56 

GaskeU C. M., Ferland G. J., 1984, PASP, 96, 393 
Georgakakis, A., Clements, D. L., Bendo, G., Rowan- 
Robinson, M., Nandra, K., Brotherton, M. S., 2009, MN- 
RAS, 394, 533 

Glikman E., Gregg M. D., Lacy, M., Helfand D. J., Becker 

R. H., White R. L., 2004, ApJ, 607, 60 
Glikman E., Helfand D. J., White R. L., Becker R. H., 

Gregg M. D., Lacy, M., 2007, ApJ, 667, 673 
Glikman E., 2011, PoS(Bashll)009 
Glikman E., et al. 2012, ApJ, 757, 51 



28 M. Rose, C. N. Tadhunter, J. Holt, J. Rodriguez Zaurin 



Goodrich R.W., 1989, ApJ, 342, 224 
Greene J.E., Ho L. C, 2005, ApJ, 630, 122 
Hao, C. N., Xia, X. Y., Mao, S., Wu, H., Deng, Z. G., 2005, 
ApJ, 625, 78 

Hazard C., Mackey A. J., Shimmins A. J., 1963, Nature, 
197, 1037 

Heckman T. M., Kauffmann G., Brinchmann J., Chariot 
S., Tremonti C., White S. D. M., 2004, ApJ, 613, 109 

Holt J., Tadhunter C., Morganti R., Bellamy M., Gonzalez 
Delgado R. M., Tzioumis A., Inskip K. J., 2006, MNRAS, 
370, 1633 

Holt J., Tadhunter C. N., Gonzalez Delgado R. M., Inskip 
K. J., Rodriguez Zaurin J., Emonts B. H. C., Morganti 
R., Wills K. A., 2007, MNRAS, 381, 611 

Holt J., Tadhunter C. N., Morganti R., 2008, MNRAS, 387, 
639 

Hopkins, P. F., Hernquist, L., Cox, T. J., Di Matteo, T., 
Martini, P., Robertson, B., Springel, V., 2005, ApJ, 630, 
705 

Hopkins P. F., Hernquist L., Cox T. J., Di Matteo T., 
Robertson B., Springel V., 2006, ApJS, 163, 1 

Howarth I. D., 1983, MNRAS, 203, 301 

Howarth I. D., Murray J., Mills D., Berry D. S., 
2004, Starlink User Note 50.24 (Swindon: PPARC), 



http://www.starlink.rl.ac.uk/star/docs/sun50.htx/sun50.htmll 



Hutchings J. P., Maddox, N., Cutri, R. M., Nelson, B.O., 

2003, AJ, 126, 63 
Jin C, Ward M., Done C, Gelbord J., 2012, MNRAS, 420, 

1825 

Jin C, Ward M., Done C, 2012, MNRAS, 422, 3268 
Jin C, Ward M., Done C, 2012, MNRAS, 425, 907 
Kewley L. J., Groves, B., Kauffmann, G., Heckman, T., 

2006, MNRAS, 372, 961 
Kim D. -C, Sanders, D. B., 1998, ApJS, 119, 
Kovacevic J., Popovic, L.C., Dimitrijevic M.S., 2010, 

MSAIS, 15, 176 
Kuraszkiewicz J., Wilkes B. J., Schmidt G., Ghosh H.,; 

Smith P. S., Cutri R., Hines D., Huff E. M., McDowell J. 

C, Nelson B., 2009, ApJ, 692, 1143 
Kuraszkiewicz J., Wilkes B. J., Schmidt G., Smith P. S., 

Cutri R., Czerny B., 2009, ApJ, 692, 1180 
LaMassa S, M., Heckman T, M., Ptak A., Martins L., Wild 

v., Sonnentrucker P., 2010, ApJ, 720, 786 
Low F.J, Cutri R.M., Huchra J.p., Kleinmann S. G., 1988, 

ApJ, 327, 41 

Marble A. R., Hines . C, Schmidt G. D., Smith P. S., 

Surace J. A., Armus L., Cutri R. M., Nelson B. O., 2003, 

ApJ, 590, 707 
Mathur S., 2000, MNRAS, 314, 17 
McLure R.J., Dunlop J. S., 2001, MNRAS, 327, 199 
Narayanan, D., Cox, T. J., Robertson, B., Dave, R., 

Di Matteo, T., Hernquist, L., Hopkins, P., Kulesa, C, 

Walker, C. K., 2006, ApJ, 642, 107 
Nelson C. H., Whittle M., 1995, ApJS, 99, 67 
Osterbrock D.E., Pogge R., 1985, ApJ, 297, 166 
Peacock J. A., 1983, MNRAS, 202, 615 
Pounds, K. A., Done, C, Osborne, J. P. 1995, MNRAS, 

277, 5 

Reyes, R., et al. 2008, AJ, 136, 2373 

Robinson T. G., Tadhunter C. N., Axon D. J., Robinson 

A., 2000, MNRAS, 317, 922 
Rodriguez Zaurin J., Tadhunter C. N., Gonzalez Delgado 



R. M. 2009, MNRAS, 400, 1139 

Rodriguez Zaurin J., Tadhunter C. N., Gonzalez Delgado 
R. M. 2010, MNRAS, 403, 1317 

Rose M., Tadhunter C. N., Holt, J., Ramos Almeida, C, 
Littlefair, S. P., 2011, MNRAS, 414, 3360 

Sarzi M., Falcon-Barroso J., Davies R. L., Bacon R., Bu- 
reau M., Cappellari M., de Zeeuw P. T., Emsellem E., 
Fathi K., Krajnovic D., Kuntschner H., McDermid R. M., 
Peletier R. F., 2006, MNRAS, 366, 1151 

Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 
525 

Schmidt M., Green R. F., 1983, ApJ, 269, 352 

Seaton M. J., 1979, MNRAS, 187, 73 

Shields G. A., Bonning E. W. 2008, ApJ, 682, 758 

Shi Y., Ogle P., Rieke G. H., Antonucci R., Hines D. C, 

Smith P. S., Low F. J., Bouwman J., Willmer C, 2007, 

ApJ, 669, 841 

Smith P. S., Schmidt G. D., Hines D. C, Cutri R. M., 

Nelson B. O., 2002, ApJ, 569, 23 
Spergel D. N., et al., 2003, ApJS, 148, 175 
Tadhunter C, Wills K., Morganti R., Oosterloo T., Dickson 

R., 2001, MNRAS, 327, 227 
Tadhunter C, Robinson T. G., Gonzalez Delgado R. M., 

Wills K., Morganti R., 2005, MNRAS, 356, 480 
Urrutia T., Lacy M., Becker R. H., 2008, ApJ, 674, 80 
Urrutia T., Becker R. H., White R. L., Glikman E., Lacy 

M., Hodge J., Gregg M. D., 2009, ApJ, 698, 1095 
Urrutia, T., Lacy, M., Spoon, H., Glikman, E., Petric, A., 

Schulz, B., 2012, ApJ, 757, 125 
Veron-Cetty M.-P., Veron P., Gongalves A. C, 2001, A&A, 

372, 730 

Webster R.L., Francis P.J., Peterson B.A., Drinkwater 

M.J., Masci F.J., 1995, Nature, 375, 469 
Weymann R. J., Morris S. L., Fohz C. B., Hewett P. C, 

1991, ApJ, 373, 23 
Whittle M., 1985, MNRAS, 213, 1 

Wilkes B. J., Schmidt G. D., Cutri R. M., Ghosh H., Hines 
D. C, Nelson B., Smith P. S., 2002, ApJ, 564, 65 

York D. G., Adelman J., Anderson J. E., Jr., Anderson S. 
F., Annis J., Bahcall N. A., Bakken J. A., Barkhouser R., 
Bastian S., Berman E., and 135 coauthors, 2000, AJ, 120, 
1579 

Zakamska N. L. et al., 2003, AJ, 126, 2144 
Zhou H., Wang T., Yuan W., Lu H., Dong X., Wang J., Lu 
Y., 2006, ApJS, 166, 128 



APPENDIX A: INDIVIDUAL PROPERTIES 

J0221 + 13 is classified as an intermediate type AGN (type 

( 2002 1 and Kuraszkiewicz et al. 



in both Smith et al. 



( |2009a,b[ l, however our classification is consistent with the 
findings of Hutchings et al. (20031: a type 2 AGN. Such a 



difference can be explained by the difference in S/N between 

the spectra used to identify the AGN type. 

J0306-05: while we find no evidence of broad line emission 



in J0306-05, Hutchings et al. (20031 classify this object as a 
type 1 AGN. 

J0312+07 is a type 1 Seyfert galaxy with J-Ks = 1.98. 
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Figure Al. The Gaussian fit leS^tKS'Mgll] AA2798,2802 absorp- 
tion and emission features. Note that both the Mgll] absorption 
and emission features overlap, making it difficult to obtain accu- 
rate measurements of the absorption/emission line properties. 



A companion object is detected on the 2D spectrum at a 
distance of ~400±80 kpc from the nucleus of J0312-)-07. 
Interestingly, the continuum of this companion appears red- 
der than the continuum of J0312-I-07. It is possible that the 
companion is interacting with the host galaxy of J0312-I-07, 
implying that J0312+07 may be a young AGN triggered in a 
recent merger. By measuring the difference in redshift in the 
Ha emission between J0312+07 and the companion galaxy, 
we have found that the companion is shifted by -240±20 
km relative to the rest-frame of J0312+07. Therefore, 
we can reasonably argue that J0312-I-07 is in an interacting 
system. 

J0400+05 is a type 1 quasar with J-Ks = 2.00. J0400-f05 
is the only object in the 2MASS sample with a detectable 



absorption feature for Mgll] AA2796,2802. Figure Al shows 



the fit to the Mgll] absorption feature. From this fit we de- 
duce that the absorption feature is blueshifted by 12Q0±160 
km s^^ relative to the rest-frame of the 0400-1-05 narrow 
line emission, and has a velocity width of 1510±120 km s^^ 
(FWHM). These properties faU short of the LoBAL prop- 



erties found for the objects in the Urrutia et al. ( 2009 1 and 



Glikman at al. ( 2012 ) samples because BAL objects are de- 
fined to by blueshifted by ^3000-25000 km s~^ and have 
FWHM > 2000 km s^ ( [Weymann et al.|l991 l. 

As well as absorption lines, the continuum of this ob- 
ject has an unusual shape which includes an unidentified, 
broad emission feature between 4200 and 5200 A in the rest- 
frame (see figure [TJ. In addition, the BLR Balmer emission 
in J0400-I-05 is significantly blue shifted by -2870±160 km 
s~^ relative to the NLR emission. Such shifts in the BLR are 
predicted for the recoils ('kicks') that result from the merger 
of two black holes ( Shields fc Bonning|2008 1 . The magnitude 
of the kick can be large: |Dain, Lousto fc Zlochowc r ( 2008| 
have computed kicks of up to 3300 km s~ , but the most 
extreme kicks require black hole binaries of equal mass and 
spin, where the spins are extremely high. We further note 
that J0400-I-05 has the highest virial black hole mass calcu- 
lated for the 2MASS sample. We note that [Hutchings et al.| 
( |2003[ ) does not detect the presence of the shifted broad line 
emission in J0400-I-05; they classify this object as a type 2 
AGN. 



ence of broad line emission in J1006-I-41, whereas we detect 
broad line emission in Ha and Fell. 

J1040+59 is a type 1 Seyfert galaxy with J-Ks = 3.02 — 
far higher than the median near-IR colour of the sample. In 
addition to its unusually red J-Ks colour, the optical con- 
tinuum of this object rises more steeply to the red than any 
other object in the sample. Its shape is reminiscent of the 
object PKS 1549-79 which also has extremely red colours 
(J-Ks=2.82; see [Bellamy et al.|2003| |Hoh et al.|2'006) . PKS 
1549-79 is believed to be a quasar which is enshrouded in a 
cocoon of natal gas and dust which is supplied by a recent 

J1040+59 



merger (Tadhunter et al. 2001 Holt et al. 



also features in the Glikman et al. (2007 



2006 1 
sample. Consis- 



tent with our results, while Glikman et al. (2007j) found that 
J1040-I-59, and many of their objects, are sufciently extin- 
guished such that the broad Balmer emission lines are not 
present at optical wavelengths, broad Paschen lines are de- 
tected in the near-IR. 
J1057-13: 



Hutchings et al. (20031 classify J1057-13 as a 



LINER/HII type o bject. However, based on the BPT dia- 
grams (see Figures 9|11 1, and presence of FHILs, we find 
clear evidence that the emission-line properties of J1057-13 
are consistent with those of a Seyfert type AGN. 
J1131 + 16 is a type 2 quasar with J-Ks = 2.15. This object 
has an unusually rich spectrum of FHILs, which are strong 
compared to the [OIII]A5007 line, and has been studied in 
detail by Rose et al.| ( |2011[ ). We believe that the unusual 
strengths of the FHILs in this object are due to a specific 
viewing angle of the far wall of the torus, coupled with a 
lack of dust on larger scales that might otherwise obscure 
our view of the torus. The spectrum of the hot dust emission 
(~1500 K) is expected to peak at wavelengths similar to 
the K-Band. The enhanced flux in this band could be the 
explanation for the red J-Ks colour observed for this object. 
J1158-30: we note that Hutchings et al. (20031 classify 



J1158-30 as a type 1 AGN, but we see no evidence of broad 
line emission in our spectra. 

J1307+23 is a type 2 Seyfert galaxy/LINER with J-Ks = 
3.31 (the reddest J-Ks colour in the sample) and has one of 
the largest redshifts in the sample (z = 0.274). It is possible 
that this object is a genuine type 2 AGN in which the AGN 
is subject to extremely high levels of extinction because the 
torus is viewed close to edge on, but due to the larger scale 
of the NLR, shows very little evidence for dust extinction 
in the narrow Balmer decrements. However, both [Smith et| 
( 2002 I and Hutchings et al. ( 2003 1 report the presence 



al 



of broad line emission in J1307-I-23, inconsistent with our 
analysis of its spectrum. On the other hand, we note that. 



like us, Kuraszkiewicz et al. (2009a,bl classify this object as 



a type 2 object. 

J1321 + 13: like J1057-13, [Hutchings et al.[ ([2003[) classify 



J1321-I-13 as a LINER/HII ty pe ob ject. Again, based on 
the BPT diagrams (see Figures 9|11 1, and the detection of 
FHILs, we find that the emission-line properties of J1321-I-13 
are consistent with those of a Seyfert type AGN. 
J1637+25 is classified as an intermediate type AGN in both 
Smith et al.|2002[ (type 1.9) and [Kuraszkiewicz et al.|2009a,b 



(type 1.5), however our classification is consistent with the 



J1006+41: Hutchings et al. ( 2003 1 does not detect the pres- 



finding of [Hutchings et al. (20031 that this object is a type 
2 AGN. 

J2124-17 is a type 1 quasar with J-Ks = 2.37. Based on 
the visual inspection of the optical spectrum, this object 
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appears to have an optical SED that is typical of UV/optical 
selected type 1 quasars. Another notable property of this 
object is that it has the strongest Fell emission of all of the 
red 2MASS sample. However, a single Gaussian fit to the 

emission line gives a velocity width FWHM of 2040±30 
km s~^ — which is higher than expected for a NLSl (FWHM 
< 2000 km s"^). 

Another interesting property of this object is that it 



has been previously identified as an ULIRG ( Kim & Sanders 



1998). Indeed, this object, along with a further 25 ULIRGs, 



is studied in more detail in [Rodriguez Zaurm, Tadhunter fc] 



Gonzalez Delgado (2009, 20101. 



A note on the comparison of classifications with 
other work. Most of the observed differences between AGN 
classifications noted in this section can be explained a by dif- 
ferences in S/N between the studies, or by differences in the 
method used to fit the Ha+[NII] blend of the observations, 
particularly for objects with weak broad lines, or classified 
as LINER/intermediate types. 



